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INTRODUCTION 


1.  SOLID  STATE  DEVICE  RESEARCH 

Depletion-mode  p+-AlInAs/InP  junction  FETs  have  been  fabricated  in  semi-insulating  InP 
using  a  selective  molecular-beam-epitaxy  technique.  Current-voltage  measurements  on  devices  with  a 
4-^m  gate  length  show  a  transconductance  of  41  mS/mm  at  zero  gate  bias,  and  RF  measurements 
indicate  a  maximum  frequency  of  oscillation  of  3.2  GHz. 

A  dual-polarization  technique  for  reducing  intermodular  on  distortion  in  integrated-optic  inter¬ 
ferometric  modulators  has  been  developed.  For  a  given  depth  of  optical  modulation,  the  intermodu¬ 
lation  distortion  can  be  significantly  reduced  by  adjusting  the  relative  amounts  of  TE  and  TM  optical 
power. 

Concurrent  fabrication  of  InGaAsP/InP  buried-heterostructure  laser  stripes  and  mirrors  has 
been  demonstrated  by  the  use  of  dry  etching  and  mass  transport.  In  comparison  to  previous  pro¬ 
cesses  in  which  the  laser  stripes  and  mirrors  are  fabricated  separately,  this  process  offers  the  advan¬ 
tages  of  simplicity,  built-in  alignment,  and  improved  control  of  device  dimensions. 

Coherent  arrays  of  multiple-stripe  InGaAs/InP  mass-transported  buried-heterostructure  lasers 
have  been  made  by  using  Y-junction  coupling.  Sharply  defined  Y-junctions  have  been  formed  by 
ion-beam-assisted  etching  to  facilitate  low-loss  coupling. 

Adsorption  on  impregnated  carbon  has  been  studied  as  a  method  to  safely  and  effectively  con¬ 
trol  arsine  and  phosphine  emissions  from  OMVPE  systems.  The  capacities  for  these  hydrides,  with 
multiple  adsorption/ oxidation  cycles,  have  been  determined  under  typical  laboratory-scale  OMVPE 
operating  conditions. 


2.  QUANTUM  ELECTRONICS 

A  reduction  in  efficiency  of  Q-switched  Nd3+  lasers  due  to  a  nonzero  lower-state  lifetime  has 
been  predicted  by  a  rate-equation  model.  For  the  1.064-^im  Nd:YAG  transition,  the  efficiency  reduc¬ 
tion  is,  at  worst,  32  percent. 

Intracavity  frequency  modulation  of  a  Nd:YAG  laser  at  very  nearly  the  cavity  mode  spacing  has 
produced  output  having  either  a  single  carrier  frequency  or  two  carrier  frequencies  corresponding  to 
adjacent  cavity  modes.  By  using  the  same  modulation  drive  frequency  for  two  Nd:YAG  lasers  (1.32 
and  1.06  ^m)  and  adjusting  the  relative  phases  of  the  drive  voltages,  the  degree  of  modulation  of 
sum-frequency  radiation  (0.589  /um)  may  be  varied  from  zero  to  twice  that  of  the  Nd:YAG  lasers. 

A  Nd:YAG  laser  system  has  been  modified  to  produce  I50-/as  macropulses,  and  frequency¬ 
doubling  experiments  have  been  carried  out  with  several  doubling  crystals  using  these  macro-pulses 
at  a  10-Hz  PRF.  In  the  best  case,  second  harmonic  generation  efficiency  was  found  to  be  limited  to 
~30  percent,  with  the  limitation  imposed  by  the  fundamental  laser  beam  configuration. 


The  time  dependence  of  amplified  spontaneous  emission  has  been  studied  using  a  simple  model. 
Experimental  measurements  of  the  gain  in  a  Ti:Al203  amplifier  affected  by  amplified  spontaneous 
emission  are  in  excellent  agreement  with  the  model. 

A  new  configuration  Faraday  isolator  for  use  with  the  TiiA^Oj  laser  has  been  constructed  using 
a  45°  birefringent  rotator  and  a  45°  Faraday  rotator.  Isolation  >27  dB  and  transmission  in  the 
forward  direction  better  than  55  percent  have  been  achieved  over  a  wavelength  range  from  700 
to  900  nm. 

3.  MATERIALS  RESEARCH 

The  infrared  absorption  of  Ir  and  IrSi  thin  films  on  Si  substrates  has  been  determined  by  trans¬ 
mission  and  reflection  measurements  over  the  wavelength  range  from  2.5  to  25  Mm-  The  existence  at 
the  Ir/ Si  or  IrSi/Si  interface  of  a  thin  boundary  layer  with  lower  absorption  than  Ir  or  IrSi  is 
indicated  by  the  thickness  dependence  of  the  absorption  and  has  been  confirmed  by  the  detection  of 
oxygen  at  the  interface  by  Auger  analysis. 

A  GaAs /  AlGaAs  light-emitting  diode  (LED)  has  been  monolithically  integrated  with  a  Si  driver 
circuit  composed  of  10  MOSFETs.  By  applying  a  stream  of  voltage  pulses  to  the  input  of  the  driver 
circuit,  the  LED  output  has  been  modulated  at  rates  exceeding  100  MHz. 

4.  MICROELECTRONICS 

An  analog  CCD  discrete  Fourier  transform  device  has  been  designed,  utilizing  a  new  split-gate 
input  structure.  The  device  performs  a  16-point  transform  with  <37-dB  harmonic  distortion  at 
5-MHz  data  rate  and  shows  greatly  improved  accuracy  between  measured  and  calculated  vector 
output  compared  with  the  earlier-reported  discrete  cosine  transform  chip. 

A  process  has  been  developed  to  uniformly  thin  and  package  large  CCD  imagers  for  use  with 
back  illumination  in  order  to  achieve  greater  green  and  blue  sensitivity.  The  noise  and  charge 
transfer  efficiency  are  not  materially  affected  by  the  process,  but  the  dark  current  is  elevated  to  ten 
times  that  of  an  unthinned  wafer. 

Various  low-temperature  dielectric  coatings  have  been  investigated  to  find  one  which  contains  a 
net  negative  charge  in  order  to  accumulate  and  stabilize  the  exposed  p-type  surface  of  a  back- 
illuminated  CCD  imager.  It  has  been  discovered  that  electron-beam-evaporated  aluminum  oxide 
contains  negative  charge  density  in  excess  of  101 1  electrons  per  square  centimeter,  a  value  adequate 
to  improve  the  response  of  the  imager  in  the  visible  and  UV,  whereas  sputtered  or  evaporated  silicon 
dioxide  or  sputtered  aluminum  oxide  has  positive  charge. 


xiv 


5.  ANALOG  DEVICE  TECHNOLOGY 

A  10-MHz  resonator  utilizing  superconductive  NbN  thin-film  electrodes  and  ceramic  thick-film 
dielectrics  has  been  built  in  order  to  demonstrate  the  feasibility  of  integrating  the  two  technologies  in 
analog  signal  processing  devices.  The  NbN  remained  superconducting  after  the  800° C  firing  of  the 
ceramic,  and  the  conductor  and  dielectric  losses  were  small  enough  to  give  a  resonator  Q  of  80. 

Calculations  of  the  phase  noise  of  10-GHz  oscillators  stabilized  with  high-temperature  super¬ 
conducting  stripline  resonators  have  been  made  and  indicate  excellent  performance.  The  residual 
noise  of  niobium-based  resonators  has  been  measured  and  exhibits  the  expected  behavior  for  offset 
frequencies  >10  Hz. 
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1.  SOLID  STATE  DEVICE  RESEARCH 

1.1  p+-Al!nAs/InP  DEPLETION-MODE  JFET 

We  previously  reported  a  selective  molecular-beam-epitaxy  (MBE)  technique  for  fabricating 
enhancement-mode  p+-AHnAs/InP  field-effect  transistors  (FETs).'  As  shown  in  Figure  1-1,  the 
basic  device  structure,  for  both  enhancement-  and  depletion-mode  FETs,  consists  of  an  ion- 
implanted  n-type  InP  channel  layer  with  a  p+-gate  layer  formed  by  selective  MBE  growth  of  Be- 
doped  AlInAs.  Details  of  the  fabrication  procedures  are  described  in  Reference  1.  In  order  to 
realize  depletion  mode  (normally  on)  FETs,  the  n-type  channel  doping  was  increased  from 
2  X  1016  cm*3  to  8  X  1016  cm-3,  and  a  2000-A  undoped  AlInAs  buffer  layer  was  incorporated  as  a 
Be  diffusion  barrier.  Indiffusion  of  Be  during  MBE  growth2  was  thought  to  contribute  to  the 
extremely  low  zero-gate-bias  drain  saturation  current  observed  in  devices  fabricated  without  a 
buffer  layer.1 


Figure  /-/.  Schematic  diagram  of  a  p+-AI/nAs/ InP  JFET  showing  the  Si* -implanted  n-InP 
and  the  MBE-grown  AlInAs  crystalline  and  polycrystalline  regions. 


For  the  present  depletion-mode  device,  Figure  1-2  shows  typical  drain  current-voltage  (I-V) 
characteristics  of  a  4-jum-gate-length  device  having  a  source-to-drain  spacing  of  10  ^m.  The 
200-/im-wide  device  exhibited  a  zero-gate-bias  saturation  current  of  40  mA  at  a  drain  bias  of  4  to 
5  V  and  had  a  pinch-off  voltage  of  about  -8  V.  The  transconductance  per  unit  gate  width  of  sev¬ 
eral  devices  across  the  wafer  was  in  the  32  to  42  mS/mm  range.  At  a  positive  gate  bias  of  1  V 
(enhancement-mode  operation),  the  saturation  current  and  transconductance,  as  shown  in 


I 


Figure  1-2.  Typical  l-V  characteristics  of  a  4  X  200-pm 
p* -AllnAsj InP  JFET  (with  buffer  layer).  The  zero-gate- 
bias  curve  is  second  from  the  lop. 


Figure  1-2,  increased  to  50  mA  and  50  mS/mm,  respectively.  Considering  the  long  gate  length 
and  large  source-drain  spacing  used,  these  values  of  scaled  transconductance  compare  very  favor¬ 
ably  with  conventionally  fabricated  InP  JFETs.3'9  Although  the  transconductance  and  pinch-off 
voltage  are  quite  reasonable,  further  improvement  could  be  expected  by  reducing  the  substantial 
gate-source  junction  leakage  current  observed  in  these  devices  (300  pA  at  pinch-off).  One  possible 
explanation  for  the  high  leakage  is  a  large  generation  rate  in  the  Alin  As  buffer  region.  Other 
possibilities  include  interfacial  trapping  centers  caused  by  a  slight  lattice  mismatch  and/or  insuffi¬ 
cient  surface-cleaning  procedures  prior  to  MBE  growth.  The  exact  cause  and  remedy  of  this 
problem  are  still  under  investigation. 

Selected  4-^m-gate-length  devices  were  cleaved  from  the  wafer  and  bonded  to  carrier  pack¬ 
ages  for  RF  characterization.  Small-signal  microwave  performance  measurements  showed  a  unity- 
power-gain  frequency  f,^  of  3.2  GHz  and  a  unity-current-gain  frequency  fT  of  2  GHz. 
Improvements  in  the  high-frequency  performance  could  be  achieved  by  scaling  down  the  gate 
length  and  employing  additional  source-drain  implants  to  reduce  parasitic  resistance. 

These  results  indicate  that  the  selective  growth  method  is  a  viable  technique  for  fabricating 
high-frequency,  high-power  junction  FETs  in  the  InP-based  materials  system. 

J.D.  Woodhouse  M.J.  Manfra 

J.P.  Donnelly  RJ.  Bailey 

1.2  REDUCTION  OF  INTERMODULATION  DISTORTION  IN 

INTERFEROMETRIC  OPTICAL  MODULATORS 

In  previous  reports,  it  was  shown  theoretically10  and  experimentally11  that  the  dynamic  range 
of  analog  optical-fiber  communication  links  employing  interferometric  external  optical  modulators 
could  be  limited  by  modulator  nonlinearities.  Specifically,  two-tone  intermodulation  (IM)  distor¬ 
tion  was  investigated.  Here,  a  dual-polarization  technique  for  reducing  IM  distortion  in  interfer¬ 
ometric  modulators  is  described. 
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The  dual-polarization  optical  link  is  illustrated  in  Figure  1-3  and  consists  of  an  optical 
source,  an  interferometric  modulator,  and  an  optical  detector.  It  is  assumed  that  the  modulator 
supports  single  TE  and  TM  modes.  A  polarizer  (or  polarization-preserving  fiber)  placed  before 
the  modulator  is  used  to  adjust  the  relative  amounts  of  TE  and  TM  input  optical  power.  For 
each  polarization,  the  output  optical  power  Pout  from  the  modulator  can  be  expressed  as 


where  Vin(t)  is  the  modulator  drive  voltage,  VK  is  the  voltage  change  required  for  full  on-off 
modulation,  6  is  a  constant  phase  bias,  and  Pin  is  the  input  optical  power.  The  detector  output 
voltage  Vout(t)  varies  linearly  with  the  sum  of  the  two  polarization  components  of  the  modulator 
output  power. 


POLARIZER 


MODULATOR 


Figure  1-3.  Analog  optical  communication  link  consisting  of  constant-power  optical 
source,  polarizer,  interferometric  modulator,  and  detector. 


To  analyze  intermodulation  distortion,  it  is  assumed  that  the  modulator  drive  voltage  Vin(t) 
consists  of  two  simultaneous  sinusoidal  signals  at  frequencies  a>j  and  to2,  each  of  amplitude  Vj^. 
As  previously  shown,10  for  small  depths  of  modulation  the  output  power  spectral  amplitudes  P|  0 
for  the  fundamentals  (a>(  and  eo2),  and  P|  2  for  the  dominant  intermodulation  signals  (2a>2  -  wt 
and  2«0|  -  <o2),  are  given  by  the  expressions 

1  ^Vm 

P|  0  =  -  —  Pin  -  sin  0  +  higher-order  terms 

^  *  TT 

3 

sin  6  +  higher-order  terms 

indicating  a  dominant  linear  and  cubic  drive-voltage  dependence  at  the  fundamental  and  1M  fre¬ 
quencies,  respectively.  Two  sets  of  equations  like  those  preceding  describe  the  TE  and  TM 
responses. 
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Typically,  an  electro-optic  modulator  can  be  designed  to  provide  different  sensitivity  for  the 
two  polarization  components  such  that 

V™  =  -vVTE 

rr  1  i r 

This  difference  in  sensitivity  can  be  exploited  to  suppress  selectively  the  dominant  cubic  term  of 
the  IM  distortion.  This  is  accomplished  by  adjusting  the  two  polarization  components  of  the 
input  optical  power,  P^E  and  P™,  and  the  phase  bias,  0TE  and  0™,  for  each  polarization  such 
that 

PTE(v™)5Sm«TE=-PTM(VTE 

For  a  LiNbOj  device,  y  is  typically  about  three;  therefore  if  the  dc  phase  biases  are  set  to 

7T  _  7 r 

0tm  =  -y ,  ^te  ■  y 

then  the  ratio  of  TM  to  TE  power  would  be  ~27  for  optimum  third-order  IM-distortion  cancel¬ 
lation.  Only  the  higher-order  components  of  the  P|  2  IM  distortion  would  remain.  Since  most  of 
the  input  power  is  now  TM-polarized,  a  consequence  of  this  technique  is  a  reduction  in  sensitiv¬ 
ity  at  the  fundamental  frequency  by  approximately  a  factor  of  7,  as  compared  to  using  only  TE- 
polarized  light.  However,  significant  increases  in  system  dynamic  range  can  be  achieved,  particu¬ 
larly  for  applications  which  have  stringent  IM-suppression  requirements  and  are  not  limited  by 
electrical  drive  power. 

The  total  normalized  electrical  output  power  of  the  link  at  the  fundamental  and  IM  frequen¬ 
cies,  as  a  function  of  the  normalized  modulator  electrical  drive  power,  is  shown  in  Figure  1-4  by 
plotting  20  log  (I  PjiE  +  PkiMl/Pin)  versus  20  log  (2Vm/VEE)  for  P10  and  P!  2  where  Pin  is  the 
total  optical  input  power.  The  dashed  lines  indicate  the  results  when  all  the  input  is  TE-polarized, 
and  the  solid  lines  show  the  dual-polarization  results  for  optimum  IM  distortion  reduction  with 
7  =  3.  A  normalized  output  power  of  0  dB  corresponds  to  full  on-off  modulation  of  the  total 
optical  power.  As  can  be  seen,  a  significant  reduction  in  IM  distortion  can  be  achieved  for  a 
given  depth  of  modulation  at  the  fundamental  frequency.  As  the  depth  of  modulation  decreases, 
the  reduction  in  IM  distortion  increases.  Likewise,  if  the  IM  suppression  is  specified,  the  funda¬ 
mental  output  power  from  the  detector  can  be  substantially  increased.  As  shown  in  Figure  1-4,  a 
9-dB  increase  in  fundamental  output  power  can  be  obtained  with  the  dual-polarization  technique 
if  a  65-dB  IM  suppression  level  is  specified. 

Experimental  results  were  obtained  by  means  of  a  LiNbOj  interferometric  modulator 
designed  for  operation  at  1.3  /urn.  Two  sinusoidal  signals  of  equal  amplitude  at  frequencies  of 
50  and  52  kHz  were  combined  and  applied  to  one  set  of  electrodes  on  the  device.  A  dc  voltage 
was  applied  to  a  second  electrode  to  adjust  the  relative  TE  and  TM  phase  bias,  and  a  polarizer 
was  placed  before  the  device  input  to  adjust  the  relative  optical  input  power  between  the  two 


4 


Figure  1-4.  Plot  of  detector  electrical  output  power  level,  20  log  {]  P^j  TE+  Pk  /  ™\  /  Plti),  for  the  fundamental, 
P/ q,  and  third-order  IM  product,  Pjj,  as  a  function  of  modulator  electrical  drive  power,  20  log  (2Vm/  Vw  T£). 
Dashed  lines  indicate  single-polarization  operation.  Solid  lines  indicate  optimized  dual-polarization  operation 
with  y  =  3. 


polarizations.  Experimental  and  theoretical  values  of  IM  suppression  for  both  single-  and  dual¬ 
polarization  operation  are  shown  in  Table  1-1  for  several  values  of  electrical  drive  power.  Theo¬ 
retical  values  of  the  normalized  fundamental  output  power  of  the  link  are  also  shown.  For  a 
given  depth  of  modulation,  the  largest  reduction  in  IM  distortion  measured  was  21  dB,  from  -46 
to  -67  dB,  at  a  normalized  output  power  of  ~-20  dB. 

L.M.  Johnson 
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TABLE  1-1 

Summary  of  Experimental  and  Theoretical  Results 


IM  Suppression  (dB) 


Drive  Power 
20  log  (2Vm/VTE) 


Fundamental 
Output  Power  (dB) 
(Theory)  (Normalized) 


TE  +TM 


-41 

-37 

-30 

-21 

-16 


TE  Only 

Exp 

Theory 

IRIii 

-66 

II 

-56 

-46 

-46 

-23 

-25 

-10 

-11 

TE  +  TM 


1.3  InGaAsP/InP  BURIED-HETEROSTRUCTURE  LASERS  WITH 
CONCURRENTLY  FABRICATED  STRIPES  AND  MIRRORS 

InGaAsP/InP  buried-heterostructure  (BH)  lasers  with  noncleaved  mirrors  are  suitable  for 
monolithic  integration  and  batch  processing  and  have  received  considerable  attention  in  recent 
years.  High-performance  lasers  with  good-quality  mirrors  have  been  demonstrated.12'16  However, 
the  BH  stripes  and  the  mirrors  were  fabricated  separately,  in  often  complicated  processes.  We 
have  developed  a  new  fabrication  process  in  which  the  laser  stripe  and  the  mirror  are  fabricated 
concurrently.  Such  a  process  not  only  is  simpler  and  more  controllable  but  also  offers  built-in 
alignment  of  the  mirror  and  the  stripe. 

The  laser-fabrication  process  involves  the  use  of  ion-beam-assisted  etching  (IBAE)  to  simul¬ 
taneously  define  the  laser  stripes  and  mirrors,  and  mass  transport  to  bury  the  active  regions  and 
smooth  the  mirrors  at  the  same  time.  A  similar  process  was  used  recently  to  demonstrate  high- 
performance  BH  lasers  with  cleaved  mirrors.17  In  contrast  to  the  previous  process,  the  mass 
transport  is  carried  out  at  a  higher  temperature  of  ~750°C  for  0.5  to  2  h,  which  is  needed  to 
completely  smooth  the  mirrors.  Shown  in  Figure  1-5  is  a  scanning-electron-microscope  (SEM) 
photograph  of  a  laser  mesa  after  mass  transport.  Typical  mesas  have  vertical  sidewalls  and  are 
3  to  4  pm  wide  except  at  the  end  regions.  The  laser  mesa  shown  has  a  tapered  end  region  with  a 
cylindrical  mirror  designed  for  reduced  output  beam  divergence.18 

The  fabricated  lasers  have  good  operating  characteristics.  Nearly  70  percent  of  the  lasers 
which  have  one  flat  etched  and  mass-transported  mirror  (with  the  other  mirror  cleaved)  have 
threshold  currents  between  16  and  26  mA.  These  devices  are  250  nm  long  and  have  active 
regions  which  are  ~0.2  Mtn  thick  and  1.5  wide.  Other  lasers  which  have  a  420-jum-long  cavity 
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Figure  1-5.  SEM  perspective  view  of  IBAE-defined 
cylindrical-mirror  laser  after  etching  and  mass  trans¬ 
port.  The  supporting  mesa  regions  on  either  side  of 
the  laser  mesa  are  also  shown. 


and  both  mirrors  etched  and  mass-transported  (with  one  mirror  flat  and  the  other  mirror  cylin- 
drically  curved)  have  threshold  currents  as  low  as  22  mA.  The  differential  quantum  efficiencies 
just  above  threshold  are  as  high  as  24  percent  and  are  similar  for  output  from  the  etched  and  the 
cleaved  mirrors.  Output  powers  >15  mW  were  measured  from  the  uncoated  etched  mirrors.  This 
performance  is  comparable  to  that  obtained  for  lasers  with  wet-chemical-etched  and  mass- 
transported  mirrors12'13  as  well  as  for  the  best  of  the  other  InGaAsP/InP  lasers  with  dry-etched 
mirrors.14"16 

The  far-field  patterns  in  the  junction  plane  of  three  devices  with  etched  and  mass-transported 
mirrors  are  shown  in  Figure  1-6.  One  device  has  a  flat  mirror  and  the  other  two  have  cylindrical 
mirrors.  The  beam  divergence  for  the  flat  etched  mirror  is  similar  to  that  for  a  cleaved  mirror. 
There  is,  however,  a  ±20  percent  device-to-device  variation  in  the  measured  beam  divergences  for 
this  wafer,  probably  due  to  variations  in  the  active-layer  thickness  and  width.  As  expected, 
cylindrical-mirror  lasers  with  a  larger  value  of  L,  the  length  of  the  tapered  region,  have  a  larger 
emitting  spot  and  thus  a  smaller  beam  divergence.  The  full-width-at-half-maximum  (FWHM) 
beam  divergence  is  only  13°  for  the  device  with  L  =  18  *tm.  The  cylindrical-mirror  lasers  can  be 
designed  to  yield  diffraction-limited  beam  divergences  by  properly  choosing  the  radius  of  curva¬ 
ture  R  of  the  mirror  with  respect  to  L.  An  approximate,  Gaussian-mode  analysis  can  be  used  to 
estimate  the  beam  divergence.  Based  on  such  an  analysis,  the  predicted  beam  divergences  are 
10.5°  for  a  device  with  R  =  12  and  L  =  18  jam,  and  19°  for  a  device  with  R  =  8  and  L  =  8  fim. 
The  predicted  divergences  are  slightly  lower  than  the  measured  values  shown  in  Figure  1-6,  since 
the  idealized  Gaussian  modes  have  the  smallest  beam  divergence. 

It  is  worth  noting  that  the  large-numerical-aperture  cylindrical  mirrors  achieved  here  cannot 
be  formed  by  wet-chemical  etching,  as  the  latter  results  in  sloped  sidewalls  and  severely  distorted 
mirror  shapes  because  of  the  crystallographic  orientation  dependence  of  the  etching.  The  present 
use  of  dry  etching  also  makes  possible  concurrent  fabrication  of  the  mirrors  and  the  stripes,  since 
the  narrow  stripe  mesas  have  vertical  sidewalls  and  result  in  narrow  transported  regions  (needed 
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Figure  1-6.  Far-field  patterns  of  IBAE-defined  etched-mirror  lasers. 


for  low  current  leakage)  even  at  the  high  mass-transport  temperature  used  for  smoothing  the  mir¬ 
rors.  In  conclusion,  the  present  concurrent  fabrication  technique  offers  advantages  such  as  built- 
in  mirror-stripe  alignment  and  controlled  fabrication  of  highly  curved  mirrors  as  demonstrated  by 
the  fabrication  of  lasers  with  small  and  predictable  beam  divergences. 

D.  Yap 
J.N.  Walpole 
Z.L.  Liau 


8 


96616-3 


1.4  Y-JUNCTION-COUPLED  InGaAsP/InP  BURIED- 

HETEROSTRUCTURE  LASER  ARRAYS 

In  recent  years,  much  effort  has  been  devoted  to  the  development  of  coherent  diode  laser 
arrays.  An  effective  method  for  coupling  these  lasers  involves  the  incorporation  of  Y-j unctions  to 
connect  the  lasing  stripes.19  This  coupling  scheme  favors  the  excitation  of  the  in-phase  array 
modes  rather  than  the  180°  out-of-phase  modes  favored  by  other  array  structures.  Here  we  report 
the  first  coherent  arrays  of  InGaAsP/InP  BH  lasers  that  are  coupled  by  Y-junctions.  In  contrast 
to  previously  reported  coherent  InGaAsP/InP  laser  arrays,20-22  the  present  arrays  are  strongly 
index-guided  and  have  the  potential  for  low  threshold  currents  and  stable,  single-lobe  outputs  at 
high  powers. 

A  schematic  illustration  of  the  Y-j  unction-coupled  array  structure  is  shown  in  Figure  1  -7(a). 
The  arrays  have  eight  stripes  on  one  end  and  nine  stripes  on  the  other  end.  The  mesa  structure 
was  fabricated  by  the  combined  use  of  IBAE  and  mass  transport  with  which  high-performance 
single-stripe  lasers  have  recently  been  demonstrated. 17  The  capability  of  the  dry-etching  technique 
to  form  narrow  grooves  with  vertical  sidewalls  without  regard  for  the  crystallographic  orientation 
is  crucial  to  the  fabrication  of  the  present  array  structure,  because  relatively  closely  spaced  mesa 
stripes  and  sharply  defined  Y-junctions  can  be  formed.  Smaller  stripe  spacings  result  in  improved 
array  fill  factors  and,  thus,  fewer  lobes  in  the  far-field  pattern.  More  sharply  defined  Y-junctions 
yield  lower  losses  and  stronger  coupling.  Figure  l-7(b)  shows  a  top  view  of  an  IBAE-defined  laser 
mesa  pattern.  The  stripes  are  2  /urn  wide  and  have  a  spacing  of  6  nm.  The  branch  angle  of  each 
Y-junction  is  10°.  Such  a  structure  cannot  be  defined  by  conventional  wet-chemical  etching  tech¬ 
niques  (which  are  dependent  on  the  crystallographic  direction),  since  the  mesa  sidewalls  would  be 
sloped  and  thus  the  Y-junctions  would  be  severely  blunted. 

The  far-field  pattern  from  the  eight-stripe  end  of  a  Y-junction-coupled  array  is  shown  in 
Figure  1-8.  This  array  has  a  threshold  current  of  300  mA  and  is  470  yum  long.  Since  typical 
threshold  currents  measured  for  400-/um-long  single-stripe  lasers  from  the  same  wafer  are  26  mA, 
we  estimate  that  the  increase  in  the  threshold  current  due  to  scattering  and  mode-mismatch  losses 
at  the  Y-junctions  is  25  to  30  percent.  Note  that  the  stripes  emit  primarily  in-phase  rather  than 
180°  out-of-phase.  The  presence  of  the  strong  narrow  lebes  in  the  far-field  pattern  indicates  that 
good  coherence  between  the  lasing  stripes  has  been  achieved. 

Spatially  resolved  emission  spectra  indicate  that  most  of  the  stripes  emit  in  several  longitudi¬ 
nal  modes.  At  most,  4  to  5  stripes  are  phase-locked  at  a  given  longitudinal  mode.  Also,  for  a 
given  longitudinal  mode  the  power  in  the  various  locked  stripes  is  not  the  same.  The  4°  FWHM 
divergence  of  the  central  far-field  lobe  shown  in  Figure  1-8  is  consistent  with  these  observations. 
Power  in  the  non-phase-locked  longitudinal  modes  contributes  to  the  baseline  seen  in  the  far-field 
pattern.  Improved  phase-locking  which  extends  over  more  stripes  can  probably  be  obtained  by 
the  use  of  more  uniform  wafers  and  improved  photolithographic  definition  of  the  stripes  and 
branches. 
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Figure  1-7.  Y-junction-coupled  laser  array:  (a)  schematic  illustration,  and 
(b)  top  view  of  a  portion  of  an  IbAE-defined  mesa. 
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Figure  1-8.  Far-field  pattern  in  the  junction  plane  of  a  Y-junction-coupled 
laser  array.  The  laser  was  excited  at  1.5  times  the  threshold  current. 


The  number  of  far-field  lobes  emitted  by  the  Y-junction-coupled  laser  array  can  be  decreased 
by  improving  the  fill-factor  of  the  array.  This  can  be  achieved  by  fabricating  cylindrical  mirrors 
at  the  ends  of  each  stripe.  The  cylindrical  mirrors  can  be  fabricated  by  IBAE  and  mass  transport 
in  the  same  process  steps  as  the  stripes  (see  Section  1.3). 

In  conclusion,  we  have  demonstrated  coherent  linear  arrays  of  mass-transported  InGaAsP/InP 
BH  lasers.  These  Y-junction-coupled  arrays  show  in-phase  emission  and  strong  phase  coherence. 

D.  Yap 
J.N.  Walpole 
Z.L.  Liau 

1.5  HYDRIDE  ADSORPTION  ON  IMPREGNATED  CARBON 

The  incomplete  pyrolysis  of  the  hydrides  arsine  (AsH3)  and  phosphine  (PH3)  in 
organometallic-vapor-phase-epitaxial  (OMVPE)  growth  systems  can  result  in  significant  environ¬ 
mental  release  of  these  highly  toxic  materials.  This  study  was  conducted  to  determine  if  adsorp¬ 
tion  of  these  materials  on  impregnated  carbon  is  a  safe  and  effective  way  to  control  emissions  of 
arsine  and  phosphine  from  OMVPE  systems. 

Carbons  impregnated  with  copper,  chromium,  and  silver  catalyze  hydride  oxidation.  For 
arsine  and  phosphine,  the  appropriate  reactions  are  thought  to  be:2-3 

2AsH3  +  402  -  As205  +  3H20 
2AsH3  +  302  -  As203  +  3H20 
2PH3  +  302  -  2H3P03 

where  the  reactions  as  shown  are  exothermic,  liberating  in  excess  of  200  kcal/gmol. 
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Existing  data  in  the  literature  could  not  be  used  because  they  apply  to  adsorption  of  very 
dilute  quantities  (ppm  level)  of  the  hydrides  in  air,  characteristic  of  large  industrial  ventilation 
systems.  OMVPE  systems,  on  the  other  hand,  typically  operate  with  relatively  low  flows  (10  to 
100  slpm)  of  a  strongly  reducing  carrier  gas  (H2)  containing  ~1  percent  hydride.  Scrubbing  the 
unpyrolyzed  hydrides  under  these  conditions  is  best  accomplished  by  a  two-step  process,  with  the 
first  being  primarily  physical  and  chemical  adsorption  of  the  hydride  on  the  carbon  during  the 
actual  OMVPE  run.  Then  the  carbon  bed  is  oxidized  to  convert  the  hydrides  to  arsenic  oxides 
and/or  phosphoric  acid  by  the  reactions  shown  above.  Chemisorption  and  oxidation  of  the 
hydride  are  necessary  to  prevent  desorption.24 

Our  studies  were  carried  out  with  a  test  cell  which  consisted  of  a  stainless  steel  canister 
which  held  400  g  carbon. +  The  canister,  7.5-cm  diam.  by  35  cm  long,  was  a  size  that  could  read¬ 
ily  be  incorporated  into  the  exhaust  plumbing  of  most  OMVPE  systems.  It  had  a  removable  top 
with  metal-seal  flange  for  replacing  the  carbon  charge  and  thermocouples  inserted  radially  into 
the  carbon  at  four  different  bed  heights.  These  proved  useful  for  tracking  the  adsorption  and 
oxidation  reactions  as  they  progressed  along  the  direction  of  flow.  The  exhaust  from  the  cell  was 
monitored  for  the  hydride  concentration. t 

Initial  studies  were  conducted  to  determine  how  bed  depth  and  flow  affected  the  carbon 
capacity  for  arsine.  Oxidation  was  then  studied  by  flowing  air  through  the  cell,  and  bed  tempera¬ 
tures  were  measured  for  different  air  flows.  During  the  third  phase  of  this  work,  the  test  cell  was 
attached  to  an  OMVPE  system  and  placed  in  the  exhaust  section  following  particulate  filters 
used  to  remove  arsenic  and  phosphorus.  In  this  configuration,  the  total  carbon  capacities  for 
arsine  and  phosphine  were  obtained  from  the  repeated  use  of  the  two-step,  adsorption-oxidation 
process. 

These  studies  showed  that  increasing  the  residence  time  of  the  arsine  in  the  bed,  either  by 
decreasing  the  flow  or  increasing  the  bed  depth,  increased  the  weight  percent  of  arsine  adsorbed 
on  the  bed.  In  all  runs,  no  arsine  was  detected  (<5  ppb)  at  the  output  of  the  test  cell  until  break¬ 
through  occurred,  and  then  concentrations  rapidly  increased  and  went  above  the  saturation  limit 
of  the  hydride  monitor  (5  ppm).  Desorption  of  arsine  was  found  to  occur  when  inert  gas  was 
used  to  purge  the  bed  before  oxidation. 

During  oxidation,  restricting  the  air  flow  through  the  cell  to  250  seem  or  less  kept  the  bed 
temperatures  below  180°C,  as  recommended  by  the  carbon  manufacturer.  Greater  air  flows  pro¬ 
duce  higher  temperatures  and  increase  the  risk  of  bed  ignition.  The  thermocouples  were  useful 
both  for  establishing  an  acceptable  air  flow  and  indicating  the  completion  of  the  reaction.  Oxida¬ 
tion  rejuvenated  the  bed  even  though  no  measurable  desorption  occurred. 


tThe  impreganted  carbon  used  in  all  the  experiments  reported  here  was  Calgon  Corp.,  Type  ASC 
12  X  30  mesh. 

tHydrides  were  sensed  with  the  MDA  Scientific  Co.,  Model  7100  monitor. 
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Figure  1-9  shows  the  results  of  the  multiple  run  tests  made  with  the  test  cell  on  the  OMVPE 
system.  In  these,  10  slpm  hydrogen  containing  1  percent  arsine  or  phosphine  were  passed  through 
the  400-g  bed  (18  cm  deep).  The  adsorption  run  was  terminated  when  a  thermocouple  5  cm  from 
the  cell  outlet  peaked  in  temperature  (typically  60°C).  With  this  as  an  endpoint  criterion,  rather 
than  hydride  breakthrough,  little  or  no  desorption  occurred  during  purging  and  oxidation,  and  a 
more  accurate  determination  could  be  made  of  the  total  hydride  adsorbed.  After  the  bed  was 
purged  with  nitrogen,  250  seem  air  was  passed  through  the  cell  until  oxidation  was  complete  (bed 
temperatures  back  down  to  below  30°C).  The  adsorption  step  was  then  repeated.  As  shown  in 
Figure  1-9,  the  bed  capacity  dropped  off  slowly,  and  six  runs  with  the  same  carbon  charge  were 
possible  in  both  the  arsine  and  phosphine  adsorption  experiments.  On  the  final  runs,  low  levels 
of  hydride  were  detected  in  the  exhaust  throughout  the  run.  The  total  carbon  adsorption  weight 
capacities  were  19  percent  phosphine  and  39  percent  arsine. 


o 

UJ 

m 
0 c 
o 
u> 
Q 
< 


UJ 

2 


O 

> 


20 

- 1 - 

1 

1 

- 1 - 

- , - , - 

D  PHOSPHINE 

15 

— 

♦ 

♦  ARSINE 

□ 

♦ 

10 

- 

□ 

a 

□ 

D 

5 

- 

♦ 

♦ 

1 

1 

_  \ 

1 

♦ 

1  1 

0 

1 

2 

3 

4 

5  6 

RUN  NUMBER 


Figure  1-9.  Volume  of  hydride  adsorbed  during  six  adsorption  /  oxidation  runs.  (The  fifth 
run  of  the  phosphine  adsorption  test  was  inadvertently  run  at  a  different  concentration  and 
the  volume  was  not  accurately  determined.  For  purposes  of  determining  the  total  capacity,  the 
volume  adsorbed  was  assumed  to  be  the  same  as  during  the  fourth  and  sixth  runs.) 

In  summary,  we  have  used  impregnated  carbon  to  safely  and  effectively  control  arsine  and 
phosphine  emissions  from  OMVPE  systems.  Concentration,  flow,  and  bed  depth  affect  the 
amount  of  hydride  that  can  be  adsorbed  before  breakthrough  occurs.  Imbedded  thermocouples 
are  useful  for  monitoring  hydride  adsorption  to  prevent  breakthrough  and  minimize  desorption. 
Oxidation  increases  total  adsorption  capacity,  but  it  must  be  controlled  to  prevent  overheating. 
After  oxidation,  hydrides  could  not  be  detected  when  the  cell  was  purged  or  when  it  was  opened 
for  carbon  removal.  A  full  account  of  these  experiments  is  presented  elsewhere.25*26 

P.L.  Green  ley  t 
S.C.  Palmateer 
S.H.  Groves 

tMIT  Environmental  Medical  Service,  Industrial  Hygiene  Office. 
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2.  QUANTUM  ELECTRONICS 


2.1  MODELING  OF  Q-SWITCHED,  FOUR-LEVEL  LASERS  WITH 
NONZERO  LOWER-LEVEL  LIFETIME 


A  simple  model  to  describe  Q-switched,  four-level  lasers  with  nonzero  lower-level  lifetime  has 
been  developed  and  used  to  evaluate  efficiency  in  Nd3+  lasers.  A  typical  assumption  in  modeling 
Q-switched  Nd^+  lasers  is  that  the  lower-laser-level  relaxation  rate  is  infinitely  fast,  but  this 
assumption  is  not  supported  by  experimental  measurements  of  this  lifetime,  which  have  ranged 
from  10  ns  to  1  yus.  Some  attempts  have  been  made  to  include  the  effect  of  the  relaxation  in 
models  of  Q-switched  Nd3+  lasers,1  but  the  energy  level  structure  of  the  Nd3+  ion  was  not  prop¬ 
erly  taken  into  account.  Our  model  shows  that  the  extraction  efficiency,  which  is  the  fraction  of 
the  upper  manifold  population  that  becomes  extracted  by  photons,  depends  on  the  ratio  of  the 
initial  upper-state  population  density  to  the  threshold  population  density,  the  ratio  of  the  lower- 
laser-level  lifetime  to  the  cavity  lifetime,  and  Boltzmann  occupation  factors  of  the  laser  transition. 
If  the  lower-level  lifetime  is  much  longer  than  the  cavity  lifetime,  the  model  predicts  that  the  effi¬ 
ciency  of  a  1.064-yum  Nd.YAG  laser  is  only  68  percent  of  that  expected  in  the  case  of  infinitely 
fast  lower-level  relaxation. 


The  energy  level  diagram  of  the  laser  and  ground-state  manifolds  for  the  1.06-^m  Nd:YAG 
transition  is  shown  in  Figure  2-1.  It  is  assumed  that  the  lower  laser  level  is  unoccupied  in  thermal 
equilibrium  because  it  is  sufficiently  above  the  ground  state.  It  is  also  assumed  that  the  crystal 
field  levels  within  each  manifold  are  in  quasi-thermal  equilibrium  at  all  times  and,  thus,  their 
populations  are  Boltzmann  distributed.  The  laser  transition  is  between  two  crystal  field  levels 
u  and  /  that  have  degeneracies  gu  and  g/  and  population  densities  Nu  and  Ny,  repectively.  The 
ratios  of  Nu  and  Ny  to  the  respective  total  manifold  populations  are  constants  fu  and  fy  deter¬ 
mined  by  the  Boltzmann  distribution. 


The  starting  point  of  the  model  is  the  laser  rate  equations 


and 


dNu 

dt 


c  a<t> 


(2-la) 


(2-lb) 


<*N,  /  gu  \  Ny 

— =  f/  Nu  -  —  Ny  cad,  -  —  (2-lc) 

dt  \  g/  /  r, 

Here  t  is  time,  a  is  the  stimulated  emission  cross  section,  ry  is  the  lifetime  of  the  lower  laser  level, 
d>  is  the  cavity  photon  density,  c  is  the  speed  of  light,  L  is  the  length  of  the  gain  medium,  L'  is 
the  cavity  length,  tr  is  the  cavity  round-trip  time,  and  8  is  the  fractional  cavity  loss  per  round 
trip.  In  Nd:YAG  all  energy  levels  have  a  degeneracy  of  2.  By  appropriate  substitutions,  these  rate 
equations  can  be  written  in  terms  of  dimensionless  quantities. 
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Figure  2-1.  Energy  level  diagram  of  the  1.064-fim  Nd:YAG  transition. 


The  rate  equations  have  been  solved  numerically.  It  has  been  assumed  that  N/(0)  =  0  and 
that  d>(0)  is  small  compared  to  the  peak  value  of  <f>( t).  The  ratio  r,  given  by  r  =  Nu(0)/Nt,  where 
Nt  is  the  threshold  population  density,  has  been  varied.  The  extraction  efficiency  i7cx  is  the  ratio 
of  the  total  number  of  stimulated  photons  to  the  total  initial  upper  manifold  population  given  by 
Nu(0)/fu.  The  number  of  stimulated  photons  is  calculated  by  determining  the  rate  of  stimulated 
emission  and  integrating  over  all  time. 

Figure  2-2  shows  i7cx(f,t)  for  the  1.064-/xm  transition  in  Nd:YAG  at  room  temperature  as  a 
function  of  r  and  r,  where  r  is  the  ratio  of  to  the  cavity  lifetime  rc  given  by  t r/6.  We  can  see 
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RATIO  OF  INITIAL 
TO  THRESHOLD  INVERSION  (r) 


Figure  2-2.  Extraction  efficiency  as  a  function  of  the  ratio  r  of  the  initial  population 
inversion  density  to  the  threshold  inversion  density  and  the  ratio  r  of  the  lower-laser- 
level  lifetime  to  the  cavity  lifetime. 


that  the  nonzero  lower-laser-level  lifetime  can  reduce  extraction  efficiency  substantially.  In  the 
worst  case  of  r  =  », 

TJex(r *°°)  =  r?ex(r.°)  fu/tfu  +  </>  (2-2) 

and  Tjex(r,oo)  is  reduced  by  68  percent  relative  to  i7ex(r,0)  in  Nd:YAG.  Typical  cavity  lifetimes  for 
Q-switched  Nd3+  lasers  are  10  ns  or  less.  Thus,  even  with  the  most  optimistic  estimate  of  lower- 
laser-level  lifetime  in  Nd:YAG,  t  >  I. 

In  summary,  modeling  shows  that  nonzero  lower-laser-level  lifetime  in  Nd3+  lasers  can  reduce 
extraction  efficiency  in  Q-switched  lasers.  It  may  be  possible  to  ameliorate  this  problem  by  find¬ 
ing  laser  materials  with  either  higher  lower-state  relaxation  rates  or  more  favorable  Boltzmann 
occupation  factors. 

T.Y.  Fan 

2.2  INTRACAVITY  FREQUENCY  MODULATION  OF  Nd:YAG 
LASER  RADIATION 

Intracavity  frequency  modulation  of  laser  radiation  has  attracted  attention  since  it  was  first 
demonstrated  in  1964.2*5  This  modulation  has  several  potential  applications  including 
spectroscopy6-9  and  optical  metrology. 10  Our  interest  in  frequency-modulated  lasers  stems  from 
our  desire  for  a  Nd:YAG  laser  that  operates  over  a  wide  spectral  range  and  yet  emits  radiation 
simultaneously  at  each  allowed  cavity  frequency  over  this  range.11*14 
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Intracavity  frequency  modulation  of  laser  radiation  may  be  accomplished  by  modulating  the 
voltage  on  an  intracavity  electro-optic  crystal  at  very  nearly  (but  not  exactly)  the  longitudinal¬ 
mode  frequency  interval  of  the  laser  cavity.  Under  ideal  conditions,  the  result  of  this  modulation 
produces  a  laser  radiation  field  given  by 

E  =  E0  cos  (wt  +  T  sin  a>mt)  (2-3) 

where  u>  is  the  carrier  frequency,  a»m  is  the  modulation  frequency,  and  T  is  the  depth  of  fre¬ 
quency  modulation.  The  depth  of  frequency  modulation  is  given  by 

T  =  (6/tt)  <dcI  (toc-a>m)  (2-4) 

where  5  is  the  depth  of  modulation  for  a  single  laser  round  trip  (proportional  to  the  modulating 
voltage  applied  to  the  electro-optic  crystal  )  and  toc  =  rrc/ L  is  the  laser  cavity  mode  spacing 
expressed  in  terms  of  the  cavity  length  L  and  the  velocity  of  light  c.  Equation  (2-3)  may  be  re¬ 
written  in  terms  of  Bessel  functions  as 

OO 

E  =  E0  £  Jn  (T)  cos[(w  +  najm)t]  .  (2-5) 

n  =  -  oo 

The  frequency  spectra  of  the  modulated  laser  consist  of  a  series  of  frequencies  a»n  separated 
by  the  modulation  frequency  aim,  and  with  powers  proportional  to  the  square  of  the  Bessel  func¬ 
tions  Jn(T).  Figure  2-3  shows  the  theoretical  spectral  power  distribution  of  frequency-modulated 
radiation  at  various  depths  of  modulation  both  for  a  single  carrier  frequency  and  for  two  carrier 
frequencies  corresponding  to  adjacent  cavity  modes.  The  frequency  extent  of  the  modulated  radi¬ 
ation  is  2 T  a»m. 

By  inserting  a  LiNb03  phase  modulation  crystal  inside  a  1.32-/*m  Nd:YAG  laser  cavity,  we 
have  obtained  frequency  modulation  of  the  laser  output  radiation.  Figure  2-4  illustrates  the  spec¬ 
tral  profiles  of  the  Nd:YAG  laser  at  various  levels  of  modulation  voltage  applied  to  the  LiNb03 
crystal.  These  spectral  profiles  were  measured  by  a  scanning  confocal  Fabry-Perot  spectrum  ana¬ 
lyzer  with  a  free-spectral  range  of  8  GHz  and  a  finesse  of  about  200.  The  depth  of  modulation  of 
each  spectral  profile  was  obtained  by  comparing  that  profile  with  theoretically  generated  spectral 
profiles.  At  zero  applied  modulation,  the  Nd:YAG  laser  operated  on  three  adjacent  longitudinal 
cavity  modes  in  the  TEMqo  mode  with  an  output  power  of  500  mW.  Applying  a  sufficiently 
strong  modulation  voltage  (>5  Vrms  )  greatly  improved  the  coherence  of  the  laser  output.  This  is 
illustrated  by  the  top  two  spectral  plots  in  Figure  2-4(a).  In  the  free-running  state  the  adjacent 
cavity  modes  operate  incoherently  relative  to  one  another,  while  in  the  frequency-modulated  state 
the  output  radiation  is  completely  coherent  and  consists  of  a  single  carrier  frequency  with  side¬ 
band  frequencies  on  both  sides.  Doubling  the  modulation  voltage  gave  a  spectral  profile  that 
could  not  be  fit  by  the  Bessel  function  spectrum  of  a  single  modulated  carrier  frequency.  How¬ 
ever,  a  good  theoretical  fit  was  possible  by  assuming  two  carrier  frequencies  corresponding  to 
adjacent  laser  cavity  modes  and  that  these  dual  carriers  were  incoherent  relative  to  one  another. 
At  higher  modulation  voltages,  the  spectral  profiles  could  usually  be  matched  to  theoretical  pro¬ 
files  by  assuming  either  single  or  dual  adjacent  carrier  frequencies.  In  many  cases,  the  laser  would 
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Figure  2-3.  Calculated  frequency-modulation  spectra  for  various  depths  of  modulation  T:  (a)  single  carrier 
frequency  to^  to  =  top  +  natm  and  (b)  two  carrier  frequencies  to p  and  top'  corresponding  to  adjacent  cavity 
modes;  <o„  =  top  +  ntow  and  to^'  =  to q  +  kum.  The  two  carrier  frequencies  are  assumed  to  be  incoherent  relative 
to  one  another.  For  the  purpose  of  this  figure,  the  sideband  intensities  of  nearly  the  same  frequency  from 
different  carrier  frequencies  are  added  incoherently. 
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RELATIVE  FREQUENCY 


Figure  2-4.  Measured  frequency-modulation  spectra  of  a  1.32- pm  Nd:YAG  laser  for  various  depths  of 
modulation  F.  The  laser  had  a  cavity  mode  spacing  of  150.000  MHz  and  was  modulated  at  149.800  MHz. 
(a)  Spectral  profiles  of  the  frequency-modulated  laser  when  operating  with  a  single  carrier  frequency.  The  laser 
free-running  spectral  profile  is  shown  in  the  uppermost  plot,  (b)  Spectral  profiles  of  the  frequency-modulated 
laser  when  operating  with  two  carrier  frequencies  corresponding  to  adjacent  cavity  modes. 
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operate  on  either  single  or  dual  carrier  frequencies  depending  sensitively  on  the  laser  cavity 
alignment.  Most  often,  the  laser  appeared  to  operate  on  dual  adjacent  carrier  frequencies. 

Figures  2-5  and  2-6  show  the  experimentally  measured  dependence  of  the  depth  of  modula¬ 
tion  on  both  modulation  voltage  and  frequency.  As  can  be  seen,  the  dependence  is  in  good  qual¬ 
itative  agreement  with  that  expected  from  Equation  (2-4). 


RELATIVE  MODULATION  VOLTAGE 


Figure  2-5.  Measured  depth  of  frequency  modulation  F 
as  a  function  of  the  relative  sinusoidal  drive  voltage 
applied  to  the  intracavity  LiNbOj  crystal. 
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Figure  2-6.  Measured  depth  of  frequency  modulation  T 
as  a  function  of  the  inverse  of  the  difference  between  the 
cavity  mode  spacing  (vc  =  atc/2ir)  and  the  modulation 
frequency  (vm  =  wJ2ir). 


An  interesting  demonstration  of  the  coherence  of  a  frequency-modulated  laser  was  obtained 
by  sum-frequency  mixing  the  output  of  two  frequency-modulated  NdrYAG  lasers,  one  operating 
at  1.06  and  the  other  at  1.32  pm.  Each  Nd:YAG  laser  was  modulated  with  the  same  depth  and 
frequency  but  with  a  180°  phase  difference.  Sum-frequency  mixing  the  radiation  from  these 
frequency-modulated  lasers  in  a  nonlinear  crystal  gives  an  unmodulated  sum  frequency.15 
Figure  2-7  shows  the  spectral  profiles  of  both  NdrYAG  lasers  and  of  the  sum  radiation  for  both 
180°  and  0°  phase  differences  between  the  two  NdrYAG  laser  modulations. 

T.H.  Jeys 
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Figure  2-7.  Spectral  profiles  of  frequency-modulated  1.06-  and  1.32-pm  Nd.YAG  lasers  and 
their  sum  radiation  at  0.589  pm.  The  modulation  depth  of  the  sum  radiation  could  be  varied 
from  zero  to  twice  that  of  the  Nd.  YAG  laser  by  adjusting  the  phase  difference  A</>  between  the 
two  laser  modulations.  The  single-frequency  result  at  Atj>  =  18CP  confirms  that  each  Nd:  YA  G 
laser  was  operating  with  a  single  carrier  frequency. 

2.3  MEASUREMENT  OF  DOUBLING  EFFICIENCY  USING 
150-ms  MACROPULSES  FROM  A  Nd:YAG  LASER 

A  long-pulse  Nd.YAG  laser  system  has  been  modified  during  the  last  three  months  to  pro¬ 
duce  mode-locked  pulse  envelopes  (macropulses)  with  durations  between  10  and  200  ms.  The  abil¬ 
ity  to  operate  with  mode-locked  pulse  envelopes  of  1  to  2  to  resemble  the  large  multibeam 
laser  system  has  also  been  maintained. 

Initial  frequency-doubling  experiments  with  this  newly  modified  system  were  carried  out  with 
150-ms  macropulses  having  the  temporal  envelope  shown  in  Figure  2-8.  This  pulse  length  is  of 
particular  interest  since  it  can  serve  as  a  pump  for  a  Ti:Al203  laser  emitting  150-ms  pulses  to 
simulate  a  diode  array.  At  this  pulse  length,  the  Nd:YAG  laser  produces  a  macropulse  energy  of 
almost  2  J.  However,  thermal  birefringence  when  the  laser  system  is  operating  at  10  Hz  produces 
a  depolarizing  effect  which  reduces  the  useful  energy  for  frequency  conversion  to  1.3  to  1.4  J  per 
pulse. 

Second  harmonic  generation  (SHG)  experiments  were  carried  out  to  establish  the  capabilities 
of  various  nonlinear  crystals  for  frequency  doubling  with  the  laser  operating  in  this  150-ms  con¬ 
figuration.  For  these  experiments  the  laser  beam  was  reduced  to  a  diameter  of  ~I  mm.  The 
results  obtained  with  the  various  crystals  are  summarized  below. 
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Two  UIO3  crystals  were  used.  At  the  highest  input  level  (1.3  J  at  10  Hz)  an  external  SHG 
efficiency  of  30  percent  was  obtained  with  an  antireflection-coated  25-mm  crystal,  and  an  external 
efficiency  of  20.3  percent  with  an  uncoated  15-mm  crystal.  These  values  correspond  to  internal 
efficiencies  of  31.2  and  25.2  percent,  respectively,  which  are  far  below  the  internal  SHG  efficien¬ 
cies  previously  obtained16  with  LilOj  using  100-/zs  Nd:YAG  macropulses. 

Three  antireflection-coated  KD*P  crystals  with  lengths  of  25,  40,  and  50  mm  were  used  in 
this  study.  At  1.3  J  they  had  internal  SHG  efficiencies  of  26.5,  18.3,  and  15.5  percent,  respec¬ 
tively.  These  values  represent  an  improvement  over  results  previously  obtained  with  a  100-//S 
pulse,  but  are  far  below  the  SHG  efficiencies  achieved  in  these  crystals  when  irradiated  with 
1-/lxs  pulses.17 

A  7.2-mm  uncoated  /3-BaB204  crystal  yielded  an  internal  SHG  efficiency  of  24  percent. 
However,  the  efficiency  had  not  saturated  at  the  highest  fundamental  input  level,  which  corres¬ 
ponded  to  ~90  MW/cm2. 

LA*P  is  a  relatively  new  nonlinear  crystal  which  has  been  described  as  superior  to  KD*P 
and  to  /J-BaB204  for  frequency  conversion.18  A  25-mm  LA*P  crystal  (90  percent  deuterated)  was 
obtained  and  tested,  with  disappointing  results.  A  maximum  of  5  percent  external  SHG  efficiency 
was  achieved  at  ~47  MW/cm2,  decreasing  to  4  percent  at  90  MW/cm2.  The  crystal  was  found  to 
be  subject  to  damage  easily  by  cleavage  and  relatively  slight  mechanical  pressure.  Further  studies 
with  this  sample  will  be  undertaken  to  establish  why  its  conversion  efficiency  falls  so  far  below 
expectations. 

A  flux-grown  4.2-mm  KTP  crystal  sustained  surface  damage  after  being  subjected  to 
~50  pulses  at  40  MW/cm2.  Prior  to  the  onset  of  damage,  an  internal  SHG  efficiency  of  23  per¬ 
cent  had  been  attained. 
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Figure  2-9.  Spatial  properties  of  (a)  Nd:YAG  laser  system  output  at  1.064  pm  and  (b)  the  second  harmonic 
generated  in  a  LilOj  crystal  at  0.532  pm.  The  light  areas  in  the  centers  are  generally  the  regions  of  most  intense 
radiation,  but  some  of  the  hot  spots  within  those  regions  appear  as  dark  islands.  The  lighter  regions  on  the 
perimeter  represent  areas  of  lower  intensity  than  in  the  intervening  dark  region. 

The  frequency  conversion  efficiencies  obtained  in  the  various  nonlinear  materials,  as  de¬ 
scribed  above,  are  uniformly  disappointing.  Even  in  the  best  case,  namely,  using  a  25-mm  LilOj 
crystal,  the  second  harmonic  output  pulse  is  limited  to  ~400  mJ.  The  relatively  low  conversion 
efficiencies  are  believed  to  be  due  to  the  multimode  nature  of  the  Nd:YAG  laser  beam. 

A  number  of  experiments  appear  to  support  this  hypothesis.  First,  a  beam-focusing  experi¬ 
ment  indicated  that  the  laser  output  is  approximately  five  times  diffraction  limited.  Next,  a  study 
of  the  beam  after  near-collimation,  reduction  to  ~l-mm  diam.  and  passage  through  apertures 
varying  from  0.4-  to  1-mm  diam.  indicated  the  presence  of  hot  spots  within  the  beam.  Finally, 
both  the  laser  beam  output  and  the  frequency-doubled  beam  were  studied  directly  and  photo¬ 
graphed,  with  the  results  shown  in  Figures  2-9(a)  and  (b),  respectively.  The  fundamental  beam 
shows  considerable  asymmetry  and  the  presence  of  hot  spots;  however,  the  most  notable  feature 
when  comparing  the  figures  is  the  presence  of  two  distinct  maxima  in  the  doubled  beam  in 
marked  contrast  to  the  fundamental.  Although  Figure  2-9(b)  was  obtained  with  a  UIO3  crystal, 
similar  results  were  obtained  with  /3-BaB204  and  KD*P.  Second  harmonic  output  energy  could 
be  transferred  from  the  upper  to  the  lower  peak  or  vice  versa  by  rotating  the  crystal,  but  a  single 


maximum  matching  that  of  the  fundamental  could  not  be  obtained.  This  indicates  the  presence  of 
diverging  phase  fronts  within  the  laser  beam  and  illustrates  the  difficulty  in  achieving  phase¬ 
matching  of  the  entire  laser  output  as  the  system  is  presently  configured. 

N.  Menyuk 
M.E.  Maclnnis 


2.4  AMPLIFIED  SPONTANEOUS  EMISSION  IN  A  Ti:AI203  AMPLIFIER 


In  recent  experimental  work,19  the  time  dependence  of  gain  in  a  TkA^Oj  amplifier  was 
measured.  It  was  shown  that  the  decay  of  gain  was  much  more  rapid  than  expected  if  only  spon¬ 
taneous  decay  was  responsible  for  depletion  of  the  upper  state.  Although  it  was  conjectured  that 
amplified  spontaneous  emission  (ASE)  might  be  responsible  for  this  difference,  the  apparent 
difficulty20-26  of  calculating  ASE  inhibited  pursuit  of  this  conjecture.  The  new  rate  equation  given 
below  shows  unequivocally  that  the  fast  decay  is  a  result  of  ASE. 


The  rate  equation  accounts  for  processes  that  directly  affect  the  total  population  inversion  in 
the  amplifier.  The  population  inversion  density  is  equal  to  the  gain  per  unit  length,  g,  divided  by 
the  stimulated  emission  cross  section  a.  The  total  population  inversion  per  unit  area  is  propor¬ 
tional  to  7  =  Sq  gdx,  where  L  is  the  length  of  the  amplifier.  The  gain  of  the  amplifier  is  e?, 
where  the  time  dependence  of  7  satisfies  the  equation 
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The  three  terms  in  this  equation  describe  the  spontaneous  emission  rate,  the  pump  rate,  and  the 
stimulated  emission  rate  caused  by  the  input  to  the  amplifier  and  by  ASE,  respectively.  Here,  the 
ASE  is  modeled  by  an  additional  input  intensity,  (HI An)  1^,  assuming  that  the  population  inver¬ 
sion  is  spatially  uniform.27  Iabs  is  the  pump  intensity  absorbed  in  the  gain  medium,  fl  is  the 
effective  solid  angle  for  ASE,  r  is  the  lifetime  of  the  excited  state,  I;n  is  the  intensity  input  to  the 
laser  amplifier,  1^,  is  the  saturation  intensity,  and  Ij^  is  the  saturation  intensity  for  the  pump. 

Equation  (2-6)  is  spatially  integrated  to  describe  the  output  of  the  amplifier.  The  inversion 
(proportional  to  7)  decays  by  spontaneous  emission  at  a  rate  of  1  /  t,  is  replenished  by  pumping, 
and  is  depleted  by  extraction  of  the  energy  either  by  a  probe  beam  or  by  ASE.  An  assumption 
that  has  been  made  in  obtaining  Equation  (2-6)  is  that  the  intensity  comes  to  a  steady  state 
rapidly  compared  with  the  population  inversion.  The  time  scale  for  the  intensity  to  reach  steady 
state  is  determined  by  the  time  required  to  traverse  the  amplifier.  As  long  as  the  net  inversion 
does  not  change  significantly  on  this  time  scale.  Equation  (2-6)  can  be  derived  from  the  usual  set 
of  spatially  dependent  Equations  (2-la,b,c).  This  assumption  is  no  longer  valid  for  high-power 
amplification  of  picosecond  and  femtosecond  pulses. 


Equation  (2-6)  also  assumes  a  two-level  model  for  the  active  medium.  The  lower  state  of  the 
gain  transition  relaxes  to  the  ground  state,  and  the  upper  state  of  the  pump  transition  relaxes  to 
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the  upper  state  of  the  gain  transition  in  times  that  are  short  compared  with  the  processes  consi¬ 
dered  above.  Such  a  model  is  applicable  to  many  important  laser  transitions.  The  rate  equation 
ignores  the  spatial  distribution  of  the  population  inversion  and  ignores  any  change  in  optical 
intensity  that  is  not  related  to  a  change  in  the  population  inversion,  for  instance,  an  absorption 
by  the  gain  medium  at  the  amplifying  wavelength.  These  assumptions  are  valid  for  the  Ti:Al2C>3 
amplifier  experiment  of  Wall  et  al.19 

In  their  experiment,  a  CW  He-Ne  laser  or  single-frequency  Ti:Al203  laser28  propagated 
almost  collinearly  with  the  pump  laser  (a  frequency-doubled  Nd:YAG  laser).  The  amplified  beam 
passed  through  a  narrowband  filter  centered  at  the  amplifying  wavelength  and  was  monitored  by 
a  photodiode.  The  photodiode  signal  was  monitored  as  a  function  of  time.  The  10-ns  pump  pulse 
is  short  compared  with  the  monitoring  time  of  4  jus  and  the  spontaneous  emission  lifetime  of 
3.15  jus.  The  probe  beam  has  an  intensity  that  does  not  significantly  deplete  the  upper  state. 

Thus,  Iabs  and  Iin  in  Equation  (2-6)  can  be  neglected  after  10  ns.  For  this  situation,  the  differen¬ 
tial  equation  is  separable  (as  long  as  O  is  not  a  function  of  time).  Thus,  the  decay  of  y  depends 
only  on  its  value  and  not  on  initial  conditions.  In  the  limit  of  no  ASE,  the  analytic  solution  is  an 
exponential,  but  in  the  general  case,  the  equation  must  be  solved  numerically. 

Figure  2-10  shows  the  measured  amplified  output  as  a  function  of  time  and  the  prediction  of 
the  rate  equation  for  a  case  with  ggL  =  y(t=0)  =  5.3;  the  agreement  is  excellent.  Without  ASE,  the 
signal  would  drop  off  half  as  fast.  The  50-ns  rise  in  the  experimental  data  is  caused  by  rise  time 
of  the  detector  and  oscilloscope  and  the  finite  duration  of  the  pump  pulse,  but  is  expected  to 
have  a  minimal  effect  on  the  data  after  50  ns. 


Figure  2-10.  The  decay  ( solid  line)  of  the  amplified  signal  at  791  nm 
in  a  Ti.  AljOj  crystal  from  Reference  10.  The  model  (dots)  gives  an 
excellent  fit  to  the  decay  with  a  solid  angle  parameter  of  0.04.  The 
experimentally  measured  initial  value  of  y  is  5.27. 
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Similarly  good  fits  of  experimental  data  have  been  obtained  for  a  series  of  different  initial  y 
between  0.5  and  5.  Each  fit  gives  a  different  value  for  the  solid  angle  parameter;  these  are  plotted 
on  a  log-log  scale  in  Figure  2-11.  For  the  largest  gains,  the  uncertainty  is  approximately  the  size 
of  the  boxes.  For  the  low  gains,  the  uncertainty  is  much  larger  because  of  the  much  lower  exper¬ 
imental  signals  and  the  reduced  contribution  of  ASE  to  gain  depletion.  For  the  lowest  gain  we 
estimate  an  uncertainty  of  approximately  a  factor  of  2  in  the  solid  angle  parameter.  In  all  the 
predictions,  the  fit  to  the  experimental  data  is  as  good  as  that  shown  in  Figure  2-10. 
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Figure  2-11.  Log-log  plot  of  solid  angle  parameters  obtained  from 
model  versus  amplifier  gain.  An  intensity-independent  solid  angle  is 
expected  at  low  intensities.  At  high  intensities  such  that  the  pump 
saturates  the  inversion,  the  solid  angle  is  expected  to  decrease  as  y\ 
shown  as  a  line  with  slope  -2  on  this  log-log  plot. 


We  can  estimate  the  solid  angle  parameter  from  the  experimental  geometry:  the  pump  beam 
diameter  is  3  mm,  and  the  1/e  absorption  length  is  5  mm.  This  gives  a  ratio  of  the  solid  angle 
effective  for  ASE  to  the  total  solid  angle,  CI/4tt,  equal  to  0.2  (see  References  20  and  22).  This 
estimate  should  be  compared  with  the  value  of  0.1  that  fits  the  experimental  results.  The  factor 
of  2  discrepancy  is  not  surprising  in  view  of  the  approximations  made  in  this  model.  At  lower 
intensities,  when  the  pump  intensity  is  well  below  the  saturation  intensity,  the  solid  angle  parame¬ 
ter  should  be  independent  of  the  pump  intensity,  as  observed. 

However,  for  pump  energies  larger  than  the  pump  beam  saturation  energy  density  of 
5  J/cm2  (corresponding  to  y  ^3),  the  absorption  length  is  expected  to  be  somewhat  longer  than 
the  5  mm  given  above,  which  would  result  in  a  somewhat  smaller  solid  angle.  The  variation  in 
solid  angle  parameter  at  different  gains  shown  in  Figure  2-11  is  ascribed  to  the  saturation  of  the 
Ti3+  ions  by  the  pump  laser  beam.  If  the  pump  beam  saturates  strongly,  then  the  length  of  the 
pumped  region  is  proportional  to  y.  The  solid  angle  should  then  vary  as  y '2.  This  is  displayed  as 
a  line  with  slope  -2  on  the  log-log  plot  of  Figure  2-11.  Good  agreement  is  obtained. 
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These  results  can  also  be  used  to  predict  gain  measurements  at  other  wavelengths.  Fig¬ 
ure  2-12  shows  a  comparison  of  the  prediction  (no  adjustable  parameters)  with  gain  measure¬ 
ments  at  633  nm;  the  agreement  is  very  good.  Also  shown  is  the  gain  if  there  were  no  ASE. 
Although  the  difference  between  these  two  curves  is  small,  the  data  fit  significantly  better  to  the 
rate  equation  including  ASE.  In  doing  these  calculations,  y  is  evaluated  for  800  nm,  then  multi¬ 
plied  by  the  ratio  of  the  gain  cross  sections  at  633  and  800  nm,  which  is  25,  to  obtain  the  signal 
intensity. 


Figure  2-12.  Decay  of  the  amplified  signal  at  633  nm  in  a 
Ti:Al2Oj  crystal  from  Reference  10.  The  model  gives  an 
excellent  fit  to  the  decay  with  no  adjustable  parameters.  Also 
shown  is  the  shape  of  the  decay  curve  if  amplified  spontane¬ 
ous  emission  is  not  included. 


The  rate  equation  approach  is  even  simpler  for  a  steady  state  calculation  where  the  left  side 
of  Equation  (2-6)  is  set  to  zero.  Evaluating  the  extraction  efficiency  as  a  function  of  small  signal 
gain  gives  good  agreement  with  the  detailed  calculations19  for  a  small-aspect-ratio  gain  region. 
Thus,  Equation  (2-6)  is  expected  to  describe  the  approximate  effects  of  ASE  on  the  gain  of  laser 
amplifiers,  including  the  possibility  of  time-dependent  pumping  and  a  time-dependent  input.  The 
limitation  imposed  by  not  considering  the  spatial  variation  of  the  inversion  by  ASE  is  probably 
not  important  for  typical  amplifier  chains  where  ASE  accounts  for  less  than  50  percent  of  the 
output  power. 

The  ease  with  which  the  effects  of  ASE  can  be  calculated  should  promote  better  understand¬ 
ing  of  its  role  in  amplifier  chains.  There  is  only  one  parameter  required  to  do  the  calculations 
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and  this  can  be  estimated  from  the  geometry.  The  simple  model  is  applied  successfully  to  an 
amplifier  with  a  small  aspect  ratio  (L/D  ~1.7).  Previous  steady  state  treatments  of  this  type21-23 
have  derivations  which  include  the  assumption  of  large  aspect  ratios  (pencil-type  geometry).  The 
rate  equation  presented  here  suggests  that  previous  work  can  be  applied  to  small-aspect-ratio 
amplifiers.  Further  testing  of  the  rate  equation  to  understand  its  limitations  is  being  done. 

In  conclusion,  a  simple  rate  equation  including  the  effects  of  ASE  gives  very  good  agreement 
with  time-dependent  experiments.  A  better  understanding  of  the  effects  of  ASE  on  experimental 
results  is  obtained. 

P.A.  Schulz 
R.L.  Aggarwal 
K.F.  Wall 


2.5  WAVELENGTH-INDEPENDENT  FARADAY  ISOLATOR 


The  dispersion  of  a  Faraday  material  can  be  compensated  by  the  dispersion  of  an  optically 
active  rotator  to  improve  the  wavelength  independence  by  more  than  an  order  of  magnitude.  This 
concept  was  realized  by  Johnston  and  Proffitt29  in  building  a  unidirectional  device  for  a  ring 
laser;  in  that  case  only  2°  rotation  was  required.  By  applying  this  concept  to  a  45°  rotation,  high 
isolation  (>30  dB)  over  a  broad  wavelength  region  is  achieved.  The  dispersion  compensation  is 
accomplished  because  both  the  Faraday  rotation  and  the  optical  activity  vary  as  the  inverse 
square  of  the  wavelength.  Such  a  Faraday  isolator  has  been  built  with  better  than  27-dB  isolation 
from  700  to  900  nm. 


The  rotation  of  either  a  Faraday  or  an  optically  active  rotator  with  45°  rotation  at  a  refer¬ 
ence  wavelength  Aref  is  given  approximately  by29 
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(2-7) 


By  choosing  the  appropriate  directions,  the  difference  between  the  rotations  in  the  Faraday 
and  the  optically  active  rotators  is  0°  in  this  approximation.  For  this  case,  perfect  isolation  is 
expected  for  a  Faraday  and  a  matched  optically  active  rotator  sandwiched  between  two  orthogo¬ 
nal  polarizers.  In  actual  practice,  because  of  nonuniformities  in  the  magnetic  field,  imperfect 
polarizers,  a  slight  mismatch  in  the  reference  wavelengths,  and  the  approximation  involved  in 
Equation  (2-7),  isolation  of  >45  dB  is  difficult  to  achieve  in  a  single  isolator.30  With  appropriate 
choice  of  material,  a  single,  fixed  isolator  can  be  used  over  a  factor  of  2  tuning  range  with 
>30-dB  isolation  (0.1  percent  backward  transmission). 


Another  important  parameter  of  a  Faraday  isolator  is  the  forward  transmission.  Here,  the 
approximate  result  is  an  important  constraint.  The  transmission  is  given  by 


T  = 


(2-8) 
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The  usable  wavelength  range  depends  on  how  much  intensity  loss  in  the  forward  direction  can  be 
tolerated.  For  50  percent  forward  transmission  loss,  the  isolator  can  be  used  for  wavelengths 
between  0.817  and  1.414  Aref.  For  a  reference  wavelength  of  780  nm,  this  constraint  implies  a 
wavelength  range  of  637  to  1100  nm,  covering  the  entire  near-infrared.  Obviously,  a  different  ref¬ 
erence  wavelength  would  shift  the  wavelength  range.  It  should  be  noted  that  the  approximation 
of  Equation  (2-7)  is  only  accurate  for  wavelengths  much  longer  than  the  absorption  edge  of  the 
material.  As  the  absorption  edge  is  approached,  a  dispersion-free  Faraday  isolator  could  be  built 
only  if  the  Faraday  and  the  optically  active  rotators  are  carefully  matched,  by  having  absorption 
edges  close  to  each  other.  In  an  actual  Faraday  isolator,  the  forward  transmission  is  also  limited 
by  the  large  number  of  surfaces  and  coatings  and  sometimes  by  absorption  in  the  Faraday  rota¬ 
tor  material.  In  the  following,  the  performance  of  such  a  Faraday  isolator  is  described. 

The  Faraday  isolator  consists  of  an  FR-5  glass  Faraday  rotator,  an  optically  active  quartz 
rotator,  and  two  crystal  polarizers.  The  Faraday  and  quartz  rotators  each  have  45°  rotation  near 
780  nm.  They  are  not  perfectly  matched  but  were  readily  available  in  our  laboratory.  The  two 
polarizers  are  aligned  perpendicular  to  each  other  so  there  is  no  backward  transmission.  There  is 
90  percent  forward  transmission  at  780  nm. 

The  backward  isolation  is  limited  by  the  problems  described  above  and  can  be  as  good  as 
the  45  dB  achieved  in  a  conventional  isolator.30  Figure  2-13  shows  that  our  design  has  isolation 
>27  dB  (backward  transmission  <0.2  percent)  between  700  and  900  nm.  Also  shown  in  Fig¬ 
ure  2-13  is  the  theoretical  limit  for  an  isolator  with  an  18-nm  mismatch  in  the  Faraday  and  the 


Figure  2-13.  Backward  transmission  through  Faraday  isolator. 
The  compensated  Faraday  isolator  has  a  backward  transmission 
(solid  line)  nearly  independent  of  wavelength.  The  slight  wavelength 
dependence  is  caused  by  a  mismatch  between  the  optically  active 
rotator  and  the  Faraday  rotator  of  ~)°  at  780  nm.  Such  a 
mismatch  gives  rise  to  an  expected  backward  transmission  (dashed 
line).  The  backward  transmission  of  a  conventional  isolator  is 
shown  by  the  dotted  line. 
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quartz  rotators.  An  independent  measurement  of  the  wavelength  mismatch  gave  a  value  of 
20  ±10  nm.  The  good  agreement  indicates  that  isolation  significantly  better  than  30  dB  should  be 
achievable  over  the  entire  wavelength  range. 

Figure  2-14  shows  the  forward  transmission  through  the  Faraday  isolator  as  a  function  of 
wavelength.  The  maximum  forward  transmission  was  ~93  percent  as  a  result  of  reflection  losses 
at  the  12  optical  surfaces  (4  on  each  polarizer  and  2  on  each  rotator).  Also  shown  is  the  theoreti¬ 
cal  forward  transmission  which  includes  the  wavelength  dependence  included  in  Equation  (2-8). 
The  actual  transmission  falls  off  faster  than  expected  from  the  theory.  Although  we  have  not 
investigated  this  effect,  it  is  likely  caused  by  the  coatings  on  the  optical  surfaces.  The  conven¬ 
tional  isolator  has  somewhat  better  transmission,  but  as  shown  in  Figure  2-13,  it  has  very  poor 
isolation  away  from  its  nominal  wavelength. 


Figure  2-14.  Forward  transmission  through  Faraday  isolator.  The 
observed  transmission  (solid  line)  is  smaller  than  expected  (dashed 
line).  The  discrepancy  may  be  caused  by  imperfect  coatings.  The 
expected  transmission  in  a  conventional  isolator  :s  shown  by  the 
dotted  line. 


In  conclusion,  a  Faraday  isolator  that  includes  an  optically  active  45°  rotator  and  a 
45°  Faraday  rotator  has  been  shown  to  give  high  isolation  over  a  very  wide  wavelength  range 
with  no  tuning.  Such  a  device  has  been  successfully  demonstrated  from  700  to  900  nm  with  over 
55-percent  transmission  and  over  27-dB  isolation.  Even  better  isolation  with  improved  optics  is 
expected. 

P.A.  Schulz 
L.J.  Belanger 
S.K.  McClung 
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3.  MATERIALS  RESEARCH 

3.1  INFRARED  ABSORPTION  OF  Ir  AND  IrSi  THIN  FILMS  ON  Si 

Silicide  Schottky-barrier  detectors  are  among  the  most  promising  sensors  for  large-scale 
monolithic  infrared  imager  arrays.  State-of-the-art  Schottky-barrier  imager  arrays  employ  PtSi 
detectors,1’2  which  have  a  cutoff  wavelength  of  ~6  pm.  With  the  objective  of  obtaining  longer 
cutoff  wavelengths,  we  are  working  on  the  development  of  IrSi  detectors.  We  recently  reported3 
the  reproducible  fabrication  of  such  detectors  with  a  cutoff  wavelength  of  ~10  pm. 

Silicide  detectors  incorporate  a  silicide-Si  Schottky  diode  that  is  formed  by  depositing  a 
noble  metal  film  on  a  p-type  Si  substrate,  then  annealing  to  convert  the  metal  to  the  correspond¬ 
ing  silicide  by  reaction  with  the  substrate.  Photodetection  is  accomplished  by  the  internal  phn- 
toemission  process,  which  includes  formation  of  hole-electron  pairs  by  infrared  absorption  in  the 
silicide  film,  transport  of  excited  holes  to  the  silicide-Si  interface,  and  finally,  transfer  of  holes 
into  the  Si  substrate  by  emission  over  the  energy  barrier  at  the  silicide-Si  interface.  The  photode¬ 
tector  quantum  efficiency  is  a  sensitive  function  of  silicide  film  thickness,  since  a  reduction  in 
thickness  increases  the  probability  that  carriers  elastically  scattered  from  the  silicide  surface  will 
be  emitted  into  the  Si,  but  also  decreases  the  optical  absorption  in  the  film.  The  optimal  thick¬ 
ness  depends  on  the  optical  characteristics  of  the  silicide.  In  the  case  of  PtSi,  experimental  values 
of  these  characteristics4-6  have  been  used  in  the  design  of  high-performance  detectors,1  for  which 
the  optimal  thickness  is  ~20  A.  In  order  to  obtain  such  values  for  Ir  and  IrSi,  which  have  not 
been  reported  previously,  we  have  measured  the  infrared  absorption  of  Ir  and  IrSi  films  on  Si 
substrates  as  a  function  of  film  thickness.  In  addition,  we  have  measured  the  absorption  of  Pt 
and  PtSi  films  on  Si  substrates  for  comparison  with  the  literature  values. 

Thin  films  of  Ir  ranging  in  thickness  from  5  to  200  A,  as  measured  with  a  quartz-crystal 
monitor,  were  deposited  by  e-beam  evaporation  of  99.9-percent-pure  Ir  on  p-type  Si  wafers  pol¬ 
ished  on  both  sides.  To  form  IrSi,  the  samples  were  removed  from  the  evaporator  and  annealed 
at  500°C  for  30  min  in  a  separate  vacuum  system.  A  similar  procedure  was  employed  to  prepare 
PtSi/ Si  control  samples.  On  the  basis  of  earlier  thickness  comparisons  made  by  transmission 
electron  microscopy,  it  has  been  assumed  that  each  silicide  film  is  twice  as  thick  as  the  metal  film 
from  which  it  was  formed. 

A  Perkin-Elmer  Model  983  spectrophotometer  was  used  to  measure  the  reflection  R  and 
transmission  T  of  the  film/ substrate  samples  over  the  wavelength  range  from  2.5  to  25  pm  for 
radiation  illuminating  the  substrate  at  near-normal  incidence,  as  shown  schematically  in  Fig¬ 
ure  3-1.  The  substrate  absorption  is  negligible  over  this  range  except  for  lattice  absorption  bands 
near  9  and  17  pm.  Outside  these  spectral  bands,  the  absorption  of  the  film,  Af,  is  equal  to  the 
absorption  of  the  film/ substrate  sample,  A  =  I  -  R  -  T.  Outside  the  lattice  bands,  the  measured 
values  of  R  and  T,  and,  therefore,  the  values  of  Af,  were  essentially  independent  of  wavelength. 
For  each  set  of  metal/ Si  and  silicide/ Si  samples,  the  R  and  T  data  obtained  at  4  pm  as  a  func¬ 
tion  of  film  thickness  were  analyzed  in  detail  in  order  to  compare  the  results  with  those  reported4 
for  PtSi  films  at  this  wavelength.  The  analysis,  which  uses  the  known  optical  constants  of 
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ABSORPTION  (%) 


Figure  3-1.  Configuration  for  measurement  of  reflection  R  and 
transmission  T  of  film  I  substrate  sample. 
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THICKNESS  OF  PtSi  (A) 

Figure  3-2.  Absorption  versus  thickness  for  PtSi  films  on  Si  substrates. 
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Si  to  evaluate  the  reflection  Rf  and  transmission  Tf  of  the  films,  yields  values  of  e'  and  e",  the 
real  and  imaginary  parts  of  the  complex  dielectric  constant  e  =  e'  +  it". 

In  Figure  3-2,  the  values  of  Af  at  4  ptm  obtained  for  the  present  PtSi  films  are  plotted 
against  film  thickness.  The  results  are  in  good  agreement  with  those  reported  by  Mooney.4  The 
theoretical  curve  shown  in  Figure  3-2  was  calculated  from  the  standard  expression7  relating 
absorption  to  thickness  by  using  the  values  e'  =  -108  and  t  -  103,  which  were  found  to  give  the 
best  least-squares  fit  to  the  values  of  Rf,  Tf,  and  Af  for  all  the  PtSi/Si  samples.  The  agreement 
between  theory  and  experiment  is  excellent.  Equally  good  agreement  was  obtained  for  the  Pt 
films. 

The  situation  is  more  complicated  for  the  Ir  and  IrSi  films.  The  values  of  Af  at  4  /am  for 
these  films  are  plotted  against  thickness  in  Figures  3-3(a)  and  (b),  respectively.  In  neither  case  is 
satisfactory  agreement  between  theory  and  experiment  obtained  by  using  a  single  pair  of  e'  and  e" 
values.  The  best  fits  obtained  by  this  procedure  are  represented  by  the  dashed  curves  in  Fig¬ 
ure  3-3.  For  each  set  of  films,  the  measured  values  of  Af  initially  fall  well  below  the  calculated 
curve.  To  account  for  the  discrepancy,  we  propose  that  a  boundary  layer  is  present  at  the  Ir-Si 
or  IrSi-Si  interface.  To  show  that  such  a  layer  could  account  for  the  observations,  we  have  calcu¬ 
lated  Af  versus  thickness  curves  on  the  basis  of  a  simple  model,  according  to  which  films  up  to  a 
certain  thickness  have  the  optical  properties  of  the  boundary  layer,  while  thicker  films  consist  of 
the  fixed  boundary  layer  and  a  variable  overlayer  of  Ir  or  IrSi,  as  shown  schematically  in  the 
insets  of  Figures  3-3(a)  and  (b).  The  values  of  e'  and  e"  for  the  boundary  layer  were  found  by  fit¬ 
ting  the  data  for  thicknesses  below  30  to  40  A,  and  these  values  were  then  adopted  in  obtaining 
e'  and  t”  values  for  the  overlayer  by  fitting  the  data  over  the  entire  thickness  range.  The  curves 
calculated  by  this  procedure,  which  are  shown  as  solid  lines  in  Figure  3-3,  are  in  very  good 
agreement  with  the  experimental  results.  The  «'  and  t"  values  for  the  boundary  layers  and  over¬ 
layers,  together  with  those  for  Pt  and  PtSi,  are  listed  in  Table  3-1. 

The  existence  of  the  proposed  boundary  layers  has  been  confirmed  by  Auger  analysis.  Fig¬ 
ure  3-4  shows  elemental  depth  profiles  for  both  as-deposited  and  annealed  Ir/Si  and  Pt/Si  sam¬ 
ples.  Oxygen  is  detected  at  the  Ir-Si  and  IrSi-Si  interfaces,  but  not  at  the  Pt-Si  and  PtSi-Si  inter¬ 
faces.  Since  Ir  reacts  very  strongly  with  oxygen,  it  is  not  surprising  that  residual  oxygen  in  the 
vacuum  system  is  incorporated  in  the  film  at  the  beginning  of  Ir  deposition.  This  gettered  oxygen 
would  remain  at  the  Si  interface  during  annealing,  since  Si  is  the  dominant  diffusion  species  dur¬ 
ing  the  formation  of  IrSi. 

The  thickness  of  the  silicide  films  used  in  Schottky-barrier  detectors  is  50  A  or  less.  Since 
the  data  of  Figure  3-3  show  that  the  presence  of  the  boundary  layer  significantly  reduces  the 
infrared  absorption  of  such  thin  films,  elimination  of  this  layer  would  improve  detector  quantum 
efficiency.  We  anticipate  that  such  improvement  can  be  achieved  by  carrying  out  the  Ir  deposi¬ 
tion  and  annealing  in  an  ultrahigh-vacuum  system  now  being  installed,  which  should  have  a  sub¬ 
stantially  lower  residual  oxygen  concentration  than  the  systems  used  in  the  present  experiments. 

C.  K.  Chen 
B-Y.  Tsaur 
M.  C.  Finn 
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AUGER  INTENSITY  (Arbitrary  Scale) 


TABLE  3-1 

Complex  Dielectric  Constants  at  4  pm 

Material 

e' 

e" 

Pt 

-90 

108 

PtSi 

-108 

103 

Ir 

Boundary  Layer 

-64 

76 

Overlayer 

-340 

490 

IrSi 

Boundary  Layer 

0 

16 

Overlayer 

-6 

54 

DEPTH  (Arbitrary  Scale) 


Figure  3-4.  Elemental  depth  profiles  obtained  by  Auger  analysis  for  (a)  Ir,  (b)  IrSi,  ( c )  Pi,  and 
(d)  PtSi  films  on  Si  substrates. 
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3.2  100-MHz  MODULATION  OF  GaAs/AlGaAs  LED  MONOLITHICALLY 
INTEGRATED  WITH  Si  DRIVER  CIRCUIT 

Monolithic  GaAs/Si  (MGS)  integration  is  potentially  an  advantageous  technology  for  cou¬ 
pling  Si  VLSI  circuits  with  GaAs/AlGaAs  optoelectronic  components  and  GaAs  high-speed  cir¬ 
cuits.  We  have  previously  reported8  such  integration  of  a  GaAs /  AlGaAs  LED  and  a  Si 
MOSFET,  with  LED  modulation  at  rates  up  to  27  MHz  accomplished  by  applying  voltage 
pulses  to  the  gate  of  the  MOSFET.  We  now  report  the  monolithic  integration  of  a  GaAs/ 
AlGaAs  LED  with  a  Si  driver  circuit  composed  of  10  MOSFETs.  Modulation  rates  of 
>100  MHz  have  been  achieved. 

The  diagram  for  the  LED/ driver  circuit  is  shown  in  Figure  3-5.  Figure  3-6  is  a  photomicro¬ 
graph  of  the  completed  circuit,  which  is  based  on  a  standard  Si  output  buffer  circuit  with  2-^m 
design  rule.  The  MOSFET  gate  oxide  thickness  is  260  A.  The  gate  widths  of  the  final  driver 
transistors,  Q9  and  Qt0,  are  412  and  376  fim,  respectively.  The  LED  replaces  the  buffer  circuit 
output  pad,  so  that  the  overall  area  remains  the  same.  The  LED  anode  is  connected  to  the  final 
driver  stage  by  Cr/Au  metallization,  and  the  cathode  is  connected  to  ground  through  the  heavily 
doped  GaAs/Si  heterointerface.  The  dimensions  of  the  unmetallized  opening  for  the  LED  output 
are  42  X  44  /xm. 


Figure  3-5.  Circuit  diagram  of  monolithically  integrated  LED /driver  circuit. 
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Figure  3-6.  Photomicrograph  of  LED/ driver  circuit. 


The  sequence  of  processing  steps  was  basically  the  same  as  the  one  used  for  LED/MOSFET 
fabrication.8  A  p  Si  wafer  slightly  misoriented  from  (100)  toward  (111)  was  heavily  implanted 
with  arsenic  to  form  n+  islands  where  LEDs  were  to  be  fabricated.  Around  each  island.  Si 
MOSFETs  were  fabricated  except  for  contact  openings  and  metallization,  using  a  poly-Si  gate 
process.  Successive  layers  of  SiC>2  and  Si3N4  were  deposited  over  the  entire  wafer  to  protect  the 
MOSFET  structures  during  GaAs/AlGaAs  epitaxy.  Openings  were  etched  in  the  Si3N4/SiOz  film 
to  expose  the  n+  islands. 

Molecular  beam  epitaxy  was  used  to  grow  the  following  layers  for  the  LED:  2-jim-thick 
n  GaAs,  0.5-/um-thick  n  AlojGag  7As,  0.3-^tm-thick  nominally  undoped  GaAs,  0.5-^m-thick 
p  Alo  3Gao  7As,  and  0.2-/um-thick  p+  GaAs.  The  growth  procedure  was  similar  to  the  one  de¬ 
scribed  previously,9  with  growth  temperatures  of  580°  C  for  GaAs  and  710°C  for  AlGaAs.  The 
material  grown  on  the  bare  Si  islands  was  single-crystal,  while  that  deposited  on  the  Si3N4/Si02 
was  polycrystalline. 

The  polycrystalline  GaAs/ AlGaAs  material  was  removed  by  wet  chemical  etching,  using  a 
photoresist  mask  to  protect  the  single-crystal  layers.  Plasma-enhanced  chemical  vapor  deposition 
was  used  to  deposit  SiNx  over  the  entire  wafer.  Contact  holes  for  the  MOSFETs  were  etched  in 
the  nitride-oxide  layer.  Metallization  was  performed  by  evaporating  and  patterning  A1  to  inter¬ 
connect  the  Si  devices.  The  LEDs  were  then  fabricated  in  the  single-crystal  layers.  Contact  open¬ 
ings  were  etched  for  the  LEDs,  and  Cr/Au  was  evaporated  to  form  ohmic  contacts  to  the  anode 
and  connect  the  anode  to  the  driver  transistors.  Finally,  A1  was  evaporated  on  the  back  of  the  Si 
wafer. 

Very  good  operating  characteristics  were  obtained  for  Si  test  devices  and  circuits  fabricated 
on  the  same  wafer  as  the  LED/driver  circuits.  For  enhancement-  and  depletion-mode  MOSFETs, 
the  respective  transconductance  values  were  about  100  and  120  mS/mm,  the  respective  threshold 
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voltages  were  0.85  and  -3  V,  the  subthreshold  slope  was  0.1  V  per  decade  of  current,  the  sub¬ 
threshold  leakage  current  was  <10'13  A//xm  gate  length,  and  the  breakdown  voltage  was  about 
7.5  V.  The  channel  mobility  estimated  from  transconductan.e  measurements  at  low  drain  bias 
was  about  540  cm2/V-s.  For  19-stage  ring  oscillators,  the  gate  delay  was  as  low  as  700  ps,  and 
the  power  dissipation  was  1.7  mW  per  gate.  These  results  show  that  inclusion  of  the  GaAs / 
AlGaAs  epitaxial  growth  step  in  the  present  processing  procedure  does  not  lead  to  significant 
degradation  of  the  Si  devices  and  circuits. 

For  dc  measurements  of  the  LED  light  output,  a  microscope  lens  with  numerical  aperture  of 
0.65  was  used  to  focus  the  light  on  a  1-cm-diam.  p-i-n  photodiode  detector.  Figure  3-7  shows  the 
output  power  as  a  function  of  dc  voltage  applied  to  the  input  of  the  driver  circuit.  The  LED  is 
turned  on  when  the  voltage  exceeds  2V.  With  the  LED  in  the  ON  state,  the  current  flow  is  about 
9  mA,  and  the  output  power  is  about  0.3  /iW.  To  modulate  the  LED  output,  a  stream  of  voltage 
pulses  was  applied  to  the  input  of  the  driver  circuit  with  the  drain  bias  at  5  V.  For  modulation 
measurements,  the  light  was  focused  with  the  microscope  lens  on  a  0.5-mm-diam.  avalanche  pho¬ 
todiode  detector.  Figure  3-8  shows  the  waveforms  for  the  input  voltage  and  light  output  when 
the  voltage  pulses  were  applied  at  a  rate  slightly  higher  than  100  MHz.  Still  higher  modulation 
rates  should  be  obtained  by  reducing  the  LED  series  resistance  and  the  MOSFET  gate  length. 

The  delay  time  of  the  light  output  with  respect  to  the  input  voltage  was  measured  by  applying 
voltage  pulses  with  a  small  duty  cycle.  The  delay  time,  which  results  mainly  from  Si  gate  delays, 
is  about  10  ns,  in  good  agreement  with  the  result  of  a  circuit  simulation. 

In  an  initial  LED  lifetime  measurement  performed  at  room  temperature  for  144  h,  the  driver 
input  and  drain  voltages  were  set  at  5  V.  The  current  through  the  LED  was  about  9  mA,  corres¬ 
ponding  to  a  current  density  of  about  550  A/ cm2.  The  light  output  gradually  decreased  with 
time,  probably  because  of  the  formation  of  dark-line  and  dark-spot  defects.  At  the  end  of  the 
run,  the  output  was  still  above  70  percent  of  its  initial  value.  This  is  a  very  encouraging  result,  in 
view  of  the  fact  that  the  GaAs/  AlGaAs  layers  have  a  defect  density  of  about  107  cm-2  and  stress 
of  more  than  109  dyne  cm'2.  It  is  anticipated  that  improvement  in  the  material  quality  will  lead 
to  greatly  increased  operating  lifetime. 

H.K.  Choi  G.W.  Turner 

J.P,  Mattia  B-Y.  Tsaur 
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1  2  3  4  5 

INPUT  VOLTAGE  (V) 


Figure  3-7. 


Dc  light  output  of  LED  as  a  function  of  driver  input  voltage. 


LIGHT  OUTPUT 
(2  mV/div) 


INPUT  VOLTAGE 
(5  V/div) 


5  ns/div 


Figure  3-8.  Waveforms  of  driver  input  voltage  (5  V/div )  and  LED  light 
output  (5  mV/div)  for  modulation  at  100  MHz. 
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4.  MICROELECTRONICS 


4.1  DISCRETE  FOURIER  TRANSFORM  CCD 

A  new  CCD  vector-matrix  product  device  which  is  capable  of  performing  a  16-point 
complex  discrete  Fourier  transform  (DFT)  has  been  designed,  and  wafers  have  been  fabricated. 

In  this  report,  we  describe  the  new  input  structure  implemented  on  the  DFT  device.  Test  results 
indicate  that  the  new  structure  yields  more  accurate  CCD  fixed-weight  multipliers  than  those  on 
the  discrete  cosine  transform  (DCT)  device  reported  earlier.1 

A  block  diagram  of  a  CCD  DFT  device  is  shown  in  Figure  4-1.  The  input  vector  x  is 
represented  by  a  complex  time  sequence.  The  output  complex  Fourier  transform  vector  X  is 
composed  of  the  DCT  and  discrete  sine  transform  (DST)  of  the  input  sequence.  Each  spectral 
output  of  the  DFT  chip  is  a  summation  of  32  products.  Each  product  of  this  summation  is  com¬ 
puted  by  a  CCD  fixed-weight  multiplier  which  is  a  single-channel  CCD  with  two  input  gates  as 


Figure  4-1.  CCD  discrete  Fourier  transform  device  concept. 
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shown  in  Figure  4-2.  If  a  potential  equilibration  input  method  is  used,  the  total  amount  of 
charge  injected  into  this  device  can  be  expressed  as^ 

Qt  =  Qs  +  Qb  (4-1) 

where  Qs  =  C0  (±A)  (Vs  -  Vmin/2  -  Vmax/2)  is  the  four-quadrant  signal  chaige  and 
Qb  =  CpA  (Vmx/2  -  Vmin/2  +  Vthgl  -  Vth  g2)  is  the  bias  charge.  In  the  expression  for  Qs,  the 
+  sign  is  used  for  a  multiplier  employing  a  normal  input  structure  and  corresponds  to  one  with  a 
positive  weight,  and  the  -  sign  is  used  for  a  device  employing  an  inverting  input  and  corresponds 
to  a  negative  weight.  In  the  normal  input  structure,  the  signal  is  applied  to  G2  and  Vmin  is  ap¬ 
plied  to  G,;  in  the  inverting  input  structure,  the  signal  is  applied  to  G,  and  V,,^  is  applied  to  G2. 
In  Equation  (4-1),  the  input  signal  amplitude  swing  is  limited  to  the  range  from  Vmin  to  Vmax. 
Also  C0  is  the  oxide  capacitance  per  unit  area,  A  is  the  input  gate  area,  and  Vth  g|  and  Vth  ^  are 
the  threshold  voltages  of  Gj  and  G2,  respectively. 


NEGATIVE  MULTIPLIER  POSITIVE  MULTIPLIER 


Figure  4-2.  CCD  fixed- weight  multiplier  with  normal  or  inverting  potential 
equilibration  input  structure. 


In  the  conventional  input  structure  shown  in  Figure  4-3,  different  levels  of  poly-Si  are  used 
for  gates  Gj  and  G2.  As  can  be  seen  in  Equation  (4-1),  if  there  is  any  threshold  nonuniformity 
between  the  poly- 1  and  -2  gates  and  if  there  is  any  process  variation  in  defining  the  G2  gate 
length,  a  computation  error  would  be  introduced  by  the  fixed-weight  multiplier.  A  new  split-gate 
input  structure  (Figure  4-3)  was  used  in  designing  the  DFT  chip  to  minimize  the  dependence  of 
the  device  performance  on  the  process  control.  The  first  and  third  gates  of  the  split-gate  structure 
are  functionally  equivalent  to  those  of  the  conventional  two-input-gate  structure,  while  the  second 
gate  is  biased  at  15  V  and  serves  as  part  of  the  storage  well.  In  this  new  design,  the  input  gate 
length  is  determined  only  by  the  pitch  of  the  poly-3  gates,  which  is  independent  of  etching  varia¬ 
tion.  In  addition,  both  the  signal  and  reference  voltages  are  applied  to  poly-1  gates;  therefore,  the 
device  performance  depends  on  the  uniformity  of  one  gate  material  only,  rather  than  the  unifor¬ 
mity  of  both  poly-1  and  -2  gates  simultaneously. 
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•  CONVENTIONAL  STRUCTURE 


i 

I 


]  POLY  1 


i  i  POLY  2 

l _ I 


•  NEW  SPLIT-GATE  INPUT  STRUCTURE 


Figure  4-3.  Conventional  and  new  split-gate  CCD  input  structures. 

A  photomicrograph  of  the  DFT  device  is  shown  in  Figure  4-4.  The  chip  size  is  7  X  9  mm.  A 
double-poly-Si,  double-metal  CCD  process  with  minimum  feature  size  of  5  pm  was  used  to  fabri¬ 
cate  the  device.  Table  4-1  presents  a  summary  of  the  CCD  fixed-weight  multiplier  performance, 
which  shows  that  using  split-gate  input  structure,  multiplier  weights  within  1  percent  of  their 
design  values  have  been  achieved.  The  accuracy  of  the  DFT  device  has  been  measured.  For  a 
given  input  vector,  each  element  of  the  output  vector  computed  by  the  device  is  digitized  and 
compared  with  its  expected  output  value.  Test  results  indicate  that  the  mean  and  standard  devia¬ 
tion  of  computation  errors  of  the  DFT  device  are  1 1  and  6  mV,  respectively.  The  maximum  out¬ 
put  of  the  device  is  ±1500  mV.  Therefore,  the  device  has  8-bit  accuracy.  The  above  test  results 
were  obtained  at  1-MHz  clock  rate.  Preliminary  test  results  indicate  that  at  10  MHz,  the  average 
computation  error  and  standard  deviation  increase  to  19.2  and  8.8  mV,  respectively.  At  10-MHz 
clock  rate,  harmonic  distortion  <-37  dB  has  been  produced  by  the  device.  Currently,  more  high¬ 
speed  tests  are  being  carried  out  and  results  are  being  analyzed. 

A. M.  Chiang  G.A.  Lincoln 

P.C.  Bennett  J.H.  Reinold 

B. B.  Kosicki  K.F.  Johnson 

R.W.  Mountain 
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Figure  4-4.  Photomicrograph  of  discrete 
Fourier  transform  CCD. 


TABLE  4-1 

Multiplier  Accuracy 

Designed  Weights 

Measured  Weights 

Conventional  Input 

Split-Gate  Input 

Normalized 

Weight 

Error 

(%) 

Normalized 

Weight 

Error 

(%) 

1.0000 

1.0000  ±  .0057 

0.00 

1.0000  ±  0066 

0.00 

0.9250 

0.9260  ±  .0062 

0.1 

0.91 50  ±  0051 

1.0 

0.7075 

0.7085  ±  .0047 

0.1 

0.7051  ±  .0068 

0.2 

0.3835 

0.3833  ±  .0033 

0.02 

0.3801  ±  .0016 

0.3 

-1.0000 

-1 .0247  ±  .0044 

2.5 

1 .0097  ±  .0067 

0.97 

-0.9250 

-0.9486  ±  .0060 

2.4 

0.9243  ±  .0027 

0.07 

-0.7075 

-0.7249  ±  .0038 

1.7 

0.71 23  ±  0051 

0.48 

-0.3825 

-0.3947  ±  .0023 

1.02 

0.3849  ±.001 3 

0.2 
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4.2  BACK-ILLUMINATED  THINNED  CCD  IMAGERS 

A  complete  mounting/ thinning/ wire-bonding  process  has  been  developed  to  fabricate  back- 
illuminated  CCD  imaging  devices  which  have  enhanced  response  in  the  blue  and  near-UV  spec¬ 
tral  region.  A  cross  section  of  the  finished  device  structure  is  shown  in  Figure  4-5. 


ILLUMINATION 


Figure  4-5.  Cross-sectional  view  of  thinned  CCD  device  structure. 


Key  elements  of  the  process  for  the  back-illuminated  imagers  are  shown  in  Table  4-2.  Each 
step  is  compatible  with  conventional  integrated  circuit  processes.  The  selection  of  the  initial  start¬ 
ing  epitaxial  material  affects  quality  of  the  final  surface  finish.  The  worst  back  surface  quality  is 
observed  on  thinned  devices  fabricated  on  wafers  of  p~  epitaxial  layers  on  p+  substrates  with 
resistivity  <0.01  fl-cm.  On  these  wafers,  the  slip  lines  due  to  the  strain  at  the  epitaxial-substrate 
interface  are  delineated  after  thinning.  No  slip  lines  are  seen  when  using  p~  epitaxial  layers  on 
either  n+  or  p+  substrates  with  resistivity  >0.01  fl-cm.  More  point  defects  are  developed  on  the 
back  surface  of  thinned  device  wafers  of  p_  epitaxy  on  n+  bulk  than  on  p'  epitaxy  on  p+  bulk.  In 
addition  to  having  fewer  point  defects,  p~  epitaxy  on  p+  wafers  are  also  preferable  to  p'  on  n+, 
because  trace  p+  substrate  material  left  on  the  thinned  back  surface  is  the  correct  type  to  enhance 
blue  response.  For  these  reasons,  p~  on  p+  epitaxial  wafers  have  been  chosen  for  most  future 
device  processing. 

After  the  CCD  imagers  are  fabricated  on  the  epitaxial  wafer,  an  oxide  layer  is  deposited 
over  the  front  of  the  wafer  to  improve  the  adhesion  between  the  wafer  and  the  glass  on  which  it 
will  be  mounted.  The  wafer  is  then  mounted  device-side  toward  the  glass  substrate  using  epoxy. 
We  have  improved  the  mounting  process  by  changing  from  the  original  soft,  low  temperature 
(100°C)  epoxy  to  a  harder,  higher  temperature  (300°C)  epoxy.  This  mounting  process  is  compat¬ 
ible  with  the  photoresist  and  wire-bonding  processes,  and  thinned  wafers  can  endure  the  hot- 
sulfuric-acid  cleaning  process  which  is  critical  for  proper  back  treatment. 
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TABLE  4-2 

Fabrication  Process  for  Backside-illuminated  Thinned  Imagers 

1.  Select  proper  p  /p+  epitaxial  material. 

2.  Fabricate  the  CCD  imager  wafers. 

3.  Mount  the  front  side  of  the  imager  to  a  glass  substrate. 

4.  Mechanically  remove  most  of  the  p+  bulk. 

5.  Selectively  etch  the  remaining  p+  bulk  and  expose  the  thin  expitaxial  layer. 

6.  Etch  silicon  and  oxide  to  expose  the  pads. 

7.  Saw  wafer  into  die,  mount  in  package,  and  wire  bond. 


The  wafer  is  then  thinned  by  first  lapping  and  polishing  the  majority  of  the  p+  substrate 
from  the  back.  The  remaining  p+  bulk  is  selectively  etched  away  to  the  thin  epitaxial  layer.  The 
number  of  defects  seen  on  the  etched  surface  is  close  to  that  of  the  original  epitaxial  wafer,  with 
the  addition  of  shallow  etch-pits  generated  from  point  defects  inside  the  p+  bulk.  The  depth  of 
these  pits  ranges  from  500  to  5000  A,  depending  on  the  vendor.  So  far,  we  have  not  observed  any 
performance  degradation  related  to  these  pits. 

A  plasma  etching  process  is  used  on  the  back  to  remove  all  Si  outside  the  active  CCD  array 
to  expose  the  A1  pads  for  wire-bonding.  The  photomask  used  to  perform  this  step  is  aligned  with 
respect  to  the  front  of  the  wafer  using  an  infrared  transmission  mask  aligner.  A  second  photo 
step  is  then  done  to  remove  the  Si02  on  top  of  the  A1  pads.  Figure  4-6  shows  a  whole  wafer 
processed  up  to  this  step.  Figure  4-7  shows  Au-wire  bonds  to  the  A1  pads  supported  by  the 
epoxy  underneath.  The  use  of  harder  epoxy  in  the  mounting  process  allows  A1  pads  to  be  wire- 
bonded  reliably  by  a  wedge-wedge  Al-wire  bonder. 

Several  back-illuminated  imagers  have  undergone  initial  electrical  and  optical  evaluation.  The 
noise  and  charge  transfer  efficiency  of  the  devices  appear  to  be  equally  as  good  as  before  thin¬ 
ning,  as  indicated  by  the  results  of  an  optical  imaging  test  shown  in  Figure  4-8.  A  resolution 
chart  has  been  imaged  onto  the  device  at  light  levels  which  produce  127  and  25  e/pixel  [Fig¬ 
ures  4-8(a)  and  (b),  respectively],  and  the  clarity  of  the  patterns  is  indicative  of  good  charge 
transfer  efficiency  at  these  low  charge  levels.  The  device  noise  is  ~15  e  rms  at  a  1-MHz  data  rate, 
which  is  the  same  as  that  of  a  typical  front-illuminated  device. 

The  most  noticeable  difference  in  the  devices  after  thinning  is  that  the  dark  current  is  highly 
nonuniform  and  up  to  10  times  higher  than  prior  to  thinning.  This  is  because  the  bare  etched 
surface  is  depleted  and  has  a  high  surface  generation  rate.  Variations  in  the  dark  current  are  the 
result  of  portions  of  the  p+  material  remaining  in  some  areas,  while  in  others  the  etchant  has 
reached  the  p~  epitaxial  material.  The  residual  p+  suppresses  the  depletion  and,  therefore,  the  sur¬ 
face  generation  rate.  Our  future  efforts  will  be  directed  at  back-surface  treatments  which  will 
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Figure  4-7.  SEM  photograph  of  wire  bonds  to  thinned  imaging  chip. 


(a) 


(b) 


Figure  4-8.  Resolution  chart  imaged  onto  thinned  CCD  at  low  light  levels: 
(a)  127  e/ pixel,  (b)  25  ej pixel. 


maintain  the  surface  in  a  state  of  accumulation  in  order  to  suppress  the  dark  current  and  ensure, 
high  quantum  efficiency.  Section  4.3  presents  results  on  the  suitability  of  thin  films  of  SiC>2  and 
A1203  for  back-surface  treatment. 


J.C.M.  Huang 
B.E.  Burke 
P.J.  Daniels 
A.R.  Forte 
D.C.  Harrison 


M. A.  Kaplan 
B.B.  Kosicki 
R.W.  Mountain 

N.  Usiak 


4.3  BACK  STABILIZATION  COATINGS  FOR  THINNED  CCD  IMAGERS 

The  response  of  thinned,  back-illuminated  CCD  imagers  could  be  improved  in  the  visible 
and  UV  by  introducing  an  electric  field  that  drives  photogenerated  electrons  deeper  into  the  Si, 
towards  the  wells  on  the  front  of  the  device.  One  means  of  creating  such  an  electric  field  would 
be  to  deposit  on  the  back  a  transparent  insulator  that  contains  a  negative  charge.  This  negative 
charge  would  induce  an  accumulation  of  holes  near  the  back,  and  the  consequent  band  bending 
would  drive  electrons  towards  the  front  surface. 

Aluminum  oxide  and  Vycor  films  about  50  nm  thick  have  been  deposited  on  p-  and  n-type 
Si  and  overlaid  with  A1  to  fabricate  MIS  diodes.  These  devices  were  tested  on  an  impedance 
bridge  and  the  measured  C-V  curves  analyzed.  Table  4-3  reviews  the  results  obtained  with  various 
deposited  layers.  Note  that  only  e-beam-deposited  A^Oj  yields  net  negative  charge.  A  negative 
charge  of  1.5  X  10* 1  e/cm2  would  induce  a  field  in  the  Si  of  about  1  X  105  V/cm,  which  should 
be  sufficient  to  increase  the  UV  response  of  the  CCD. 
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TABLE  4-3 

Fixed  Charge  in  Deposited  Oxides 


Oxide 

Deposition 

Method 

,A>x  _ 

1011  e/cm2 

Thermal  Oxidation 

3.1 

E-beam  (No  02) 

12 

Vycor 

(7910) 

Sputter  (No  02) 

45 

ai2o3 

Sputter  (No  02) 

32 

Al203 

E-beam  (No  02) 

-5 

Al203 

E-beam 

(4  X  10'5  Torr  02) 

-1.1 

It  is  interesting  to  note  that  the  addition  of  O2  to  the  e-beam  system  during  deposition 
decreases  the  amount  of  negative  charge.  This  is  opposite  to  what  might  be  expected  from  ther¬ 
modynamic  arguments,  and  suggests  that  the  presence  of  oxygen  alters  the  deposition  kinetics  of 
AI2O3  so  that  a  more  perfect  oxide  layer  is  deposited.  This  conclusion  is  supported  by  the  higher 
breakdown  voltage  found  for  films  deposited  in  the  presence  of  02.  In  addition  to  breakdown  at 
high  fields,  the  diodes  also  exhibit  marked  hysteresis  in  bias  stress  tests.  A  positive  bias  on  the  A1 
electrode  causes  a  negative  shift  in  the  fixed  charge,  while  a  negative  bias  causes  a  positive  shift, 
leading  eventually  to  a  net  positive  fixed  charge  in  the  AI2O3.  This  hysteresis  should  not  be  a 
cause  for  concern  in  the  devices  while  in  operation,  since  the  device  will  not  be  at  elevated 
temperature  or  under  an  applied  electric  field.  Also,  this  shift  in  net  charge  will,  whether  positive 
or  negative,  anneal  out  over  a  period  of  several  days  or  weeks  at  room  temperature,  and  the 
initial  negative  charge  will  return.  Further,  the  first  lot  of  samples,  which  has  been  unstressed, 
has  retained  negative  charge  for  over  five  months. 

The  nature  of  the  mobile  charges  in  the  A12C>3  films  is  not  clear,  as  the  defect  species  in  bulk 
AI2O3  has  not  been  definitively  established.  It  is  unlikely  that  alkali  ions  are  the  mobile  species 
since  deposited  films  of  AI2O3  are  used  to  block  Na  migration  into  Si02.  The  fixed  charge  indu¬ 
ces  a  field  in  the  oxide  layer  that  is  about  1  X  105  V/cm  and  this  corresponds  to  a  drop  across 
the  oxide  of  about  0.5  V.  Such  a  voltage  drop  does  not  cause  any  motion  of  charge  in  bias  stress 
tests  and  is  below  the  breakdown  voltage. 

It  appears,  therefore,  that  deposition  of  AI2O3  will  produce  a  stable  accumulation  of  holes 
on  the  illuminated  surface  of  a  thinned  imager  and  thereby  improve  the  response  of  the  back 
imager. 
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The  first  CCDs  (both  packaged  and  unpackaged  devices)  are  now  in  the  process  of  being 
coated  with  AI2O3.  The  thickness  chosen  is  76  nm,  so  that  the  oxide  layer  will  act  as  an  antire' 
flection  coating  and  increase  the  photon  flux  at  k  =  500  nm  by  about  50  percent. 

J.A.  Gregory 
R.W.  Mountain 
R.V.  Gray 
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5.  ANALOG  DEVICE  TECHNOLOGY 

5.1  PROCESS  INTEGRATION  OF  SUPERCONDUCTING  THIN  FILMS 
AND  DIELECTRIC  THICK  FILMS 

Compact,  superconductive  analog  signal  processing  devices  with  multigigahertz  bandwidths 
have  been  constructed  at  Lincoln  Laboratory.  One  of  the  most  important  elements  in  the  archi¬ 
tecture  of  these  devices  is  delay  lines.  The  low  RF  losses  associated  with  superconductors  allow 
electromagnetic  delay  lines,  tens  of  meters  long,  to  be  built  without  appreciable  loss,  but  their 
packing  density  is  limited  by  crosstalk  between  adjacent  lines.  This  crosstalk  is  minimized  for 
striplines,  where  it  falls  exponentially  with  the  distance  between  lines  normalized  to  the  distance  b 
between  ground  planes  (Figure  5-1).  For  crosstalk  between  40  and  60  dB,  the  required  separation 
is  between  1.5  and  2  b.  Hence,  thinner  dielectrics  are  desired  because  they  allow  a  greater  pack¬ 
ing  density  and  more  delay  per  substrate. 


Figure  5-1.  Basic  structure  for  the  stripline  using  an  integrated  superconductive-thin-film/ 
dielectric-thick-film  approach. 

The  minimum  thickness  of  the  dielectric  is  determined  by  allowable  losses.  To  maintain  the 
same  characteristic  impedance,  the  width  of  the  line  must  be  reduced  in  proportion  to  the  thick¬ 
ness  of  the  dielectric.  Under  these  circumstances  the  attenuation  of  the  lines  is  inversely  propor¬ 
tional  to  the  thickness.  This  imposes  a  lower  limit  on  the  thickness  of  the  dielectric. 

Most  of  our  current  devices  use  125-jxm-thick,  2-in-diam.  Si  wafers  for  a  maximum  delay  of 
40  ns.  This  delay  has  been  extended  to  100  ns  with  the  use  of  50-Mm-thick  wafers,  but  this 
represents  the  limit  of  the  present  technology,  since  thinner  or  larger  substrates  are  too  fragile  to 
handle  and  difficult  to  manufacture.  A  design  that  maximizes  delay,  however,  would  require 
thinner  and  larger  substrates.  For  example,  for  high-quality  Nb  films,  an  optimum  design  leads 
to  the  following  parameters:  14-Mm-wide  lines,  distance  between  ground  planes  equal  to  20  pm, 
and  distance  between  lines  equal  to  32  pm.  We  assume  a  dielectric  constant  of  10.  For  this 
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optimized  design,  a  total  delay  of  1  ps  can  be  accommodated  in  a  3-in  substrate,  with  a  maxi¬ 
mum  loss  of  1  Np  (8.69  dB)  at  10  GHz  and  a  maximum  crosstalk  of  50  dB  between  adjacent 
lines.  For  other  design  parameters,  other  dielectric  thicknesses  will  be  optimum,  but,  in  general, 
for  most  applications  a  dielectric  thickness  between  10  and  25  pm  is  desirable. 

An  alternative  method  to  fabricate  these  lines  is  an  integrated  technology  using  thin-film 
superconductors  and  thick-film  dielectrics.  In  this  approach,  shown  schematically  in  Figure  5-1,  a 
superconducting  ground  plane  is  first  deposited  on  a  support  substrate,  followed  by  a  thick-film 
dielectric.  A  second  superconductive  thin  film  is  then  deposited  on  the  dielectric  layer  and  pat¬ 
terned,  forming  the  center  conductor  of  the  delay  line.  This  is  followed  by  a  second  layer  of  the 
dielectric  and  the  second  ground  plane. 

There  are  many  difficulties  involved  with  developing  a  technology  to  fabricate  this  structure. 
The  deposition  method  for  the  dielectric  must  provide  uniform  layers,  low  dielectric  losses,  and 
be  compatible  with  the  superconductive  thin  films.  Thin-film  methods  for  deposition  of  dielectrics 
are  usually  limited  to  the  deposition  of  films  up  to  2  pm  thick  and  usually  produce  fairly  lossy 
films.  Low-loss  dielectric  bulk  ceramics  can  be  produced,  but  they  usually  require  high  firing 
temperatures  (>1300°C).  Also,  screen  printing,  the  common  method  for  formation  of  thick 
ceramic  films,  can  only  reliably  produce  films  thicker  than  50  pm.  These  films  are  usually  fired  at 
temperatures  higher  than  850°C  in  oxidizing  atmospheres,  which  would  damage  the  underlying 
superconducting  film.  Recently,  thick-film  pastes  that  can  be  fired  in  N2  and  are  capable  of  pro¬ 
ducing  thick  films  with  tan  6  ~10~3  at  room  temperature  became  available.  These  films  contain 
low-melting-point  glass  mixed  with  a  ceramic  powder  to  lower  the  firing  temperature  of  the 
ceramic. 

We  have  utilized  one  of  these  thick-film  pastes  to  study  the  feasibility  of  producing  thin-film- 
superconductor/ thick-film-dielectric  integrated  circuits.  It  was  our  hope  that  at  lower  operating 
temperatures  the  losses  associated  with  these  ceramic  thick  films  would  be  lowered  enough  to 
warrant  their  use  in  practical  analog  signal  processing  devices. 

The  basic  structure  used  as  a  test  vehicle  for  the  development  of  this  technology  is  illustrated 
in  Figure  5-2.  An  LC  resonator  is  used  to  identify  possible  problems  with  the  process  and  at  the 
same  time  evaluate  the  dielectric  loss  at  low  temperatures. 

The  superconductor  used  in  this  structure  is  NbN  because  it  is  much  less  reactive  than  Nb. 
Thin  film:,  of  NbN  are  deposited  by  reactive  sputtering  in  an  Ar-5%  N2  atmosphere  on  a 
2-in-diam.  oxidized  Si  wafer  using  an  RF  diode  system.  The  NbN  layer  used  to  define  the  bot¬ 
tom  electrode  and  the  inductor  is  5000  A  thick. 

The  thin  dielectric  layer  between  the  NbN  and  the  ceramic  layer  shown  in  Figure  5-2  has 
two  functions.  As  an  interface  layer,  it  enhances  the  adhesion  of  the  ceramic  layer  to  the  super¬ 
conductor.  As  a  diffusion  layer,  it  prevents  the  reaction  of  the  ceramic  with  the  NbN  and  the 
oxidation  of  the  latter  during  the  firing  cycle  for  the  ceramic.  Candidate  materials  for  this  layer 
were  selected  based  on  their  thermodynamic  stability.  Among  the  oxides  more  stable  than  Nb205, 
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NbN  TOP  ELECTRODE 


NbN  GROUND  PLANE 


Figure  5-2.  Exploded  view  of  the  L-C  resonator  test  structure. 


we  chose  SiO,  MgO,  and  A1203  because  of  their  known  compatibility  with  NbN.  A  series  of 
high-temperature  (annealing)  experiments  were  carried  out  to  test  the  adequacy  of  these  layers, 
leading  us  to  choose  a  multilayer  dielectric.  First,  a  thin  layer  of  Si02  (500A)  is  deposited  on  the 
NbN  without  exposing  it  to  the  air.  The  two  films  are  then  etched  utilizing  reactive  ion  etching 
(RIE).  After  patterning,  a  second  layer  of  Si02  (7000A)  and  a  layer  of  MgO  (I000A)  are  depos¬ 
ited.  This  combination  was  found  to  minimize  the  degradation  of  the  superconducting  thin  film. 

The  thick-film  paste  used  for  this  structure  is  Electro  Science  Laboratory  (ESL)  4906,  which 
is  N-j-firable  and  is  reported  by  the  manufacturer  to  have  a  tan  8  of  10'3  at  iow  frequency 
(1  kHz)  and  at  room  temperature.  Numerous  other  pastes  were  also  examined  but  were  found  to 
be  incompatible  with  some  step  of  the  fabrication  procedure.  In  order  to  form  films  of  the 
desired  thickness,  a  spin-on  process  was  developed.  The  paste  was  diluted  in  thinner  provided  by 
ESL  and  spun  on  the  patterned  wafer.  Films  10  to  20  pm  thick  could  be  easily  produced  using 
this  method.  A  liftoff  procedure  was  used  to  define  the  ceramic  area.  We  used  a  thick-film  pho¬ 
topolymer  that  was  laminated  onto  the  wafer  using  the  same  technology  as  for  metal  plating.  The 
photoresist  was  then  patterned  using  conventional  photolithography.  Dielectric  paste  was  then 
spun  onto  the  photoresist-patterned  wafer  and  dried.  The  liftoff  was  accomplished  by  dissolution 
of  the  photoresist  in  a  5  percent  solution  of  potassium  hydroxide. 

The  firing  of  the  thick-film  paste  was  carried  out  in  two  stages.  The  burnout  of  the  binder 
was  carried  out  in  air  at  450°C  for  1  h  using  a  box  muffle  furnace.  Previous  tests  had  shown 
that  the  oxide-protected  NbN  did  not  deteriorate  after  this  treatment. 

We  found  that,  even  using  extensive  N2  flushing  of  the  furnace,  N2  firing  at  850°C  would 
destroy  the  superconductivity  of  the  NbN  films.  Auger  analysis  showed  that  this  was  caused  by 
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oxidation  of  the  NbN  by  residual  oxygen  in  the  furnace.  Consequently,  a  vacuum-firing  process 
was  implemented.  For  this  vacuum-sintering  process,  the  ceramic-coated  wafers  were  placed  in  a 
quartz  tube  which  was  pumped  to  10'7  Torr  using  an  ion  pump.  The  quartz  tube  was  then  placed 
in  a  furnace  which  was  ramped  up  to  850°C  in  3  h  and  was  held  at  this  temperature  for  30  min. 
During  the  sintering,  the  outgassing  of  the  film  exceeded  the  pumping  capacity  of  the  ion  pump 
and  a  sorption  pump  was  used  for  the  last  part  of  the  sintering  sequence.  The  pressure  at  this 
stage  was  lO^  Torr.  The  furnace  w;as  then  allowed  to  cool  to  room  temperature.  A  test  made 
with  this  annealing  procedure  showed  no  significant  suppression  of  Tc  for  the  NbN  films.  (Typi¬ 
cal  Tc’s  of  15  K.  were  obtained  on  large-area  samples,  and  a  <2-percent  drop  occurred  upon 
annealing.)  Finally,  a  thick  film  of  NbN  (~2  Mm)  was  deposited  on  the  back  of  the  Si  wafer  and 
on  the  top  of  the  thick  film.  The  counterelectrodes  for  the  capacitors  were  patterned  using  RIE. 

The  device  was  packaged  and  contacts  were  made  to  the  two  lower  capacitor  electrodes  of 
the  resonator  using  ultrasonically  bonded  A1  wires.  The  packaged  device  was  tested  using  an 
impedance  analyzer  in  a  four-probe  configuration.  The  device  was  mounted  at  the  end  of  a  probe 
and  loaded  into  a  liquid-He  dewar.  Four  !-m-long  coaxial  cables  connected  the  package  to  the 
analyzer.  The  absolute  value  of  the  impedance  and  the  phase  angle  were  measured,  from  which 
the  quality  factor  (Q)  for  the  LC  circuit  was  determined. 

Two  devices  were  tested.  They  have  resonant  frequencies  of  6.1  and  5.6  MFIz  and  Qs  of  70 
and  28,  respectively.  These  values  of  Q  are  smaller  than  expected  from  the  manufacturer  data.  It 
is  not  clear  at  this  time  if  the  higher-than-expected  losses  of  the  capacitors  are  due  to  the  modi¬ 
fied  firing  schedule  or  to  undetected  damage  to  the  NbN.  In  new  structures,  measurements  of  the 
changes  of  the  Q  with  temperature  will  help  to  clarify  the  loss  mechanism. 

S.C.  Wong 
A.C.  Anderson 


5.2  HIGH-TC  SUPERCONDUCTING  RESONATORS 

The  recent  discovery  of  the  new  metal-oxide  high-transition-temperature  superconducting 
materials1  has  stimulated  interest  in  many  applications,  including  high-quality-factor  resonators 
which  could  be  operated  at  liquid  N2  temperatures  (77  K).  Such  resonators  would  be  the  basis 
for  construction  of  low-phase-noise  oscillators  for  application  to  radar  and  communication  sys¬ 
tems.  The  use  of  superconductors  is  of  interest  because  extremely  high  Qs  and  high  power  levels 
are  possible.  Qs  of  106  for  thin-film  stripline  resonators  operating  at  1  GHz  and  I05  at  10  GHz 
are  achievable  with  the  standard  materials  such  as  Nb  at  liquid  He  temperatures.  Even  higher  Qs 
have  been  achieved  using  superconducting  cavities,  but  the  stripline  structure  is  preferred  because 
it  is  compact,  easy  to  fabricate,  and  can  be  integrated  with  other  planar  circuits. 

The  resonator  Q  is  determined  by  the  RF  surface  resistance  of  the  superconducting  material 
which,  although  low,  is  not  identically  zero  as  it  is  at  dc.  The  finite  RF  resistance  is  a  result  of 
the  nonzero  mass  of  the  superconducting  electron  pairs;  an  ac  voltage  results  from  this  inertia 
which  in  turn  drives  the  lossy  normal  electrons  which  are  present  at  any  nonzero  temperature. 
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Few  measurements  of  the  RF  surface  resistance  have  been  made  in  the  new  high-Tc  materials, 
and  as  yet  there  is  no  reliable  theoretical  calculation.  As  discussed  below,  predictions  of  Q  pres¬ 
ently  must  rely  upon  measured  data  which  depend  strongly  upon  the  material  processing  and 
thus  should  not  be  considered  definitive. 


However,  if  low-loss  characteristics  can  be  achieved  reproducibly  at  gigahertz  frequencies 
with  the  new  materials,  then  devices  of  practical  importance  can  be  considered.  Detailed  calcula¬ 
tions  of  oscillator  phase  noise  will  be  presented,  including  such  factors  as  amplifier  flicker  noise 
and  power-handling  capability  of  the  superconductor.  Calculations  using  a  standard  two-fluid 
model  for  the  high-Tc  materials  based  upon  the  measured  bulk  low-frequency  properties  show 
that  a  Q  of  108  should  be  possible  at  1  GHz  at  50  K..  However,  the  validity  of  the  two-fluid 
model  has  not  been  confirmed,  and  measured  values  of  resistance  are  more  than  an  order  of 
magnitude  higher  than  those  calculated  from  standard  superconductivity  theory. 

Before  calculating  the  phase  noise  of  an  oscillator  stabilized  with  a  superconducting  resona¬ 
tor,  we  must  first  estimate  the  Q  expected  for  the  stripline  structure  and  the  maximum  power 
that  the  resonator  can  support.  In  Figure  5-3  we  show  the  cross  section  of  the  proposed  stripline 
resonator.  The  two  substrates  are  of  an  appropriate  low-loss  dielectric  material,  e.g.,  MgO,  1  mm 
thick.  The  center  conductor  is  2.54  mm  wide.  For  this  configuration,  the  conductor-limited  qual¬ 
ity  factor,  Qc,  obtainable  with  a  material  having  a  surface  resistance  Rs  in  ohms  at  a  frequency  f 
in  hertz  is 


Qc  =  3.73  X  10-9  X  f/Rs 
and  the  maximum  input  power  is 

™  z0i ; 

"  8Il(1-Il)Ql 


(5-D 


(5-2) 


where  Z0  is  the  impedance  of  the  stripline,  Ic  the  maximum  current  for  the  center  conductor,  IL 
the  insertion  loss,  and  QL  the  loaded  Q.  For  the  chosen  dimensions,  Z0  is  28  ft. 


We  assume  Rs  =  4. 1  X  10*4  ft  at  10  GHz  and  50  K.,  obtained  by  extrapolating  the  results 
from  the  Stanford/ UCLA  measurements,2  done  at  100  GHz,  to  lower  frequencies  assuming  an  f2 
dependence.  For  calculating  Ic,  we  assume  a  critical  current  of  1  X  107  A/cm2  at  50  K,  a  conser¬ 
vative  assumption  in  view  of  the  results  of  the  Stanford  group.  For  these  values,  and  assuming 
an  insertion  loss  of  15  dB,  we  obtain  a  Q  of  150,000  at  5  GHz  and  a  maximum  input  power  of 
32.4  W.  The  corresponding  values  for  a  center  frequency  of  10  GHz  are  Q  =  75,000  and 
P  =  64  W.  These  power  levels  are  maximum  input  powers  and  not  values  of  dissipated  power. 


In  the  calculation  of  Q  and  critical  current,  we  used  a  stripline  configuration  for  the  trans¬ 
mission  line.  Stripline,  besides  eliminating  radiation  losses,  can  be  dimensioned  to  maximize  the 
Q  and  input  power.  Other  structures,  notably  microstrip  and  coplanar  lines,  could  also  be  used. 
Microstrip  lines  have  the  advantage  that  they  require  only  a  single  substrate  and  hence  are  less 
sensitive  to  mechanical  vibrations.  The  radiation  losses  can  be  minimized  by  careful  design  and 
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Figure  5-3.  Stripline  resonator  structure:  (a)  cross  section;  (b)  top  view  of  center  conductor. 


packaging.  Coplanar  lines  have  the  additional  advantage  that  films  have  to  be  deposited  on  only 
one  side  of  the  substrate.  For  the  same  thin-film  conductor  properties,  coplanar  lines  are  lossier 
and  can  carry  less  current  than  a  microstrip  or  stripline,  but  they  may  be  adequate  for  our 
purposes. 


Having  calculated  the  Q  and  power,  one  can  predict  the  values  of  phase  noise  to  be  expected 
from  a  feedback  oscillator  stabilized  by  a  superconducting  resonator,  using  the  standard  Leesom 
model.3  The  stability  of  an  oscillator  to  be  used  in  a  Doppler  radar  or  communication  system  is 
best  analyzed  in  the  frequency  domain  where  one  considers  the  FM  noise  sidebands  of  the  oscil¬ 
lator.  This  quantity  is  called  the  single-sideband  phase  noise  and  is  given  by  the  following 
expression: 


L»  =  10  log 


N2 


4Q2w2> 


GFkT 

r*r 


+ 


(5-3) 


where  L(o>)  is  the  noise  power  relative  to  the  carrier  in  a  1-Hz  bandwidth  (given  in  dBc/Hz),  Pc 
is  the  power  at  the  output  of  the  amplifier,  G  is  the  loop  gain,  f  is  the  amplifier  noise  figure,  to  is 
the  offset  frequency  in  rad/s,  Q  is  the  loaded  resonator  quality  factor,  cu0  is  the  oscillator  funda¬ 
mental  frequency,  a  is  the  empirically  determined  flicker-noise  constant,  k  is  Boltzmann’s  con¬ 
stant,  T  is  the  absolute  temperature,  and  N  is  the  frequency-multiplication  factor.  The  flicker 
(1/f)  noise  can  originate  in  either  the  amplifier,  or  the  resonator,  or  both.  Since  the  flicker  noise 
is  not  well  understood  theoretically,  it  must  be  empirically  determined.  We  can  make  reasonable 
estimates  for  the  amplifier  contribution,  but  for  the  resonator  we  must  rely  on  experimentally 
determined  values.  This  is  a  very  important  point.  Our  projections  of  phase  noise  that  follow 
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assume  that  the  contribution  to  the  flicker  noise  from  the  resonator  is  small  compared  with  the 
amplifier  contribution,  but  the  phase  noise  at  offsets  of  less  than  about  5  kHz  is  dominated  by 
the  flicker  noise.  If  the  flicker  noise  contribution  of  the  superconducting  resonator  is  larger  than 
that  of  the  amplifier,  then  the  actual  phase  noise  will  be  greater  than  the  following  predictions. 
The  importance  of  the  resonator  flicker  noise  has  led  us  to  undertake  the  measurements  described 
below. 

Figure  5-4  shows  the  projected  phase  noise  at  10  GHz  for  a  superconducting-resonator- 
stabilized  oscillator  for  two  different  sets  of  parameters.  One  curve  represents  an  oscillator  oper¬ 
ating  at  5  GHz,  frequency-doubled  to  obtain  10  GHz.  The  other  parameters  are  Q  =  150,000, 

G  =  15  dB,  F  =  6  dB,  P  =  25  dBm  (0.3  W),  and  a  =  4  X  10'12.  This  value  of  a  is  that  normally 
assumed  for  a  GaAs  FET  amplifier.  The  second  curve  represents  an  oscillator  operating  at 
10  GHz  with  a  Q  -  75,000  and  the  other  parameters  unchanged  from  the  first  curve.  The  reduc¬ 
tion  in  Q,  when  raising  the  frequency  from  5  to  10  GHz,  comes  from  the  frequency  dependence 
of  the  surface  resistance  of  superconductors.  The  parameter  values  chosen  for  these  curves  are  for 
illustrative  purposes  and  should  not  be  considered  final.  We  have,  for  instance,  assumed  room- 
temperature  operation  for  the  amplifier,  but  when  it  is  operated  at  the  cryogenic  temperatures, 
noise  performance  would  be  improved.  We  have  also  assumed  an  oscillator-loop  power  of 
25  dBm  because  this  value  should  be  easily  supported  by  a  relatively  simple  and  inexpensive 
amplifier.  As  discussed  above,  the  power-handling  capacity  of  the  resonator  itself  is  considerably 
higher  than  the  25  dBm  used  in  the  calculation.  Still  better  phase-noise  performance  would  be 
obtained  with  higher  power  levels,  but  this  higher  performance  may  not  be  necessary  and  may 
not  justify  the  added  expense  of  a  higher-power  amplifier.  In  addition,  the  insertion  loss  of  the 
resonator  is  a  design  parameter  which  we  are  free  to  choose  to  optimize  overall  performance.  For 
instance,  it  may  be  useful  in  the  future  to  design  the  resonator  so  that  a  single  transistor  with  a 
10  dB  gain  provides  the  necessary  loop  gain  for  oscillation.  This  is  entirely  possible  to  do  with  a 
slight  degradation  of  loaded  Q.  Such  an  option  may  be  attractive  if  the  unloaded  Q  can  be  made 
higher  than  our  conservative  projections. 

Also  shown  in  Figure  5-4  are  the  phase-noise  curves  for  quartz  crystal  and  SAW  oscillators. 
Clearly,  the  superconducting  oscillator  provides  better  performance  than  any  other  technology 
and  is  able  to  meet  the  needs  of  real  radar  systems  with  a  comfortable  margin. 

Because,  as  discussed  above,  the  flicker  noise  of  the  superconducting  resonator  is  unknown 
and  because,  even  for  the  standard  superconducting  materials,  there  are  no  relevant  experimental 
data,  we  have  made  preliminary  measurements  of  the  residual  noise  of  an  all-Nb  superconducting 
resonator  fabricated  with  the  stripline  structure  described  above.  The  devices  have  a  loaded  Q  of 
~2.5  X  I04,  a  fundamental  resonant  frequency  of  1.15  GHz,  and  were  operated  at  4  K  with  the 
resonator  in  a  liquid-He  bath.  The  measurements  were  made  by  passing  a  constant  frequency 
tone,  generated  by  a  synthesizer,  through  the  resonator  under  test,  and  beating  the  output  of  the 
resonator  against  the  tone.  The  power  spectrum  of  this  baseband  output  is  the  resonator  contri¬ 
bution  to  the  I/f  noise,  the  last  term  in  parentheses  of  Equation  (5-3). 
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Figure  5-4.  Calculated  single-sideband  phase  noise. 

The  results  are  shown  in  Figure  5-5.  The  bandwidth  of  this  resonator  is  40  kHz,  so  one  sees 
in  the  data  the  synthesizer  phase  noise  for  offset  frequencies  >40  kHz.  For  offset  frequencies 
<1  kHz,  the  rise  in  the  data  is  the  result  of  the  resonator.  It  is  not  clear  how  much  of  the  noise 
is  due  to  fundamental  flicker  noise  in  the  resonator,  how  much  is  due  to  vibration  in  the  resona¬ 
tor  or  liquid-He  dewar  probe,  and  how  much  is  due  to  temperature  fluctuations.  The  noise  below 
10-Hz  offset  no  longer  follows  a  1/f  frequency  dependence,  so  that  it  is  probably  the  result  of 
either  thermal  or  mechanical  instabilities.  (This  resonator  design,  however,  was  made  for  material 
evaluation  and  not  for  oscillator  application,  so  it  is  not  surprising  that  performance  is  affected 
by  vibration.) 

Between  10  Hz  and  1  kHz  the  curve  closely  follows  a  1/f  dependence.  If  the  1/f  dependence 
between  10  Hz  and  1  kHz  is  typical,  the  resonator  flicker  noise  will  be  comparable  to  the  GaAs 
amplifier  contribution  and  will  not  severely  degrade  the  performance  predicted  in  Fig¬ 
ure  5-4. 


D.E.  Oates 
A.C.  Anderson 
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j  (f)  SINGLE  SIDEBAND  PHASE  NOISE 
(dBc/Hz) 


Figure  5-5.  Measured  residual  noise  for  all-Nb  superconducting  resonator  at  1.15  GHz. 
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